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ON THE BENDING OF A RAY THROUGH A 
SPHERICALLY SYMMETRIC ATMOSPHERE 

Harvey G. Safren 

Laboratory for Theoretical Studies 


ABSTRACT 

An equation is derived which explicitly gives the 
path of a ray through a spherically symmetric atmos- 
phere. The refraction of the ray is assumed to be gov- 
erned by the laws of geometric optics. The equation in- 
volves the arrival angle of the ray at the ground; this 
angle is obtained by solving an auxiliary equation. The 
latter equation is not explicitly solvable, but approximate 
solution formulas are obtained and a numerical solution 
method is described. Some of the approximate formulas 
are shown to be accurate to better than one percent down 
to elevation angles of ten degrees or less. The refrac- 
tivity profile includes both the troposphere and ionosphere, 
and is arbitrary within a very large class of functions. 
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ON THE BENDING OF A RAY THROUGH A 
SPHERICALLY SYMMETRIC ATMOSPHERE 


INTRODUCTION 

To determine the orbit of a satellite with great precision it is necessnry to 
correct ground-based radar observations for the effects of refraction the 
atmosphere. In theory it is possible to compute these corrections very accurately, 
but in practice it is not. The atmosphere is so complex, and its behavior so 
variable, that one cannot hope to have a very complete description of it at the 
time the radar observations are made. The best one can ordinarily hope to do is 
to assume a model of the atmosphere which is greatly simplified, but still real- 
istic enough so that the corrections computed from it can be regarded as 
reasonably accurate. In this report the main simplifying assumption made is 
that the atmosphere is spherically symmetric. It is also assumed that the track- 
ing frequencies used are well above 100 megahertz, so that the effects of the 
earth’s magnetic field are not too important, and that the variation of refractive 
index with height is gradual enough so that the laws of geometric optics can be 
assumed to hold. With these assumptions the problem of computing the arrival 
angle becomes manageable. (The computation of range and range-rate errors 
will be treated in a later report.) It is not necessary to assume any specific 
law for the variation of the refractive index with altitude. 

The mathematical formulation of the problem follows the method outlined 
in Ref. 1 (pp. 120-122), but the derivation given here is more complete. The 
basic idea is to make use of the differential equation for a ray path. This equa- 
tion does not give the ray path explicitly, because it uses the path length along 
the ray as the independent variable, but it is easily integrated for a spherically 
symmetric medium. The result is the formula known as Bouguer's Law, or the 
spherical form of Snell’s Law: nR sin 6 = constant, where n is the index of 
refraction at the distance R from the center of the earth, and 0 is the angle be- 
tween the ray path and the radius vector R. This relation also does not explicitly 
give the ray path in cartesian or polar coordinates, but it is easy to deduce from 
by geometrical reasoning, a first-order differential equation for the ray path 
in terms of R and the central angle </■ (see Fig. 1). This equation can be im- 
mediately integrated, but the integrated equation still contains the unknown con- 
stant in Bouguer's law. To find the constant we Impose the condition that the ray 
path must join the satellite to the ground station. This gives an equation for b , 
the elevation angle error (i.e., the angle betv n the ray path and the slant patli 
at the ground station). Unfortunately appears in a definite integral and it 
doesn’t seem possible to explicitly solve for it. ' But it is possible to deduce a 
considerable amount of information about b as a function of the altitude and 

g 

elevation angle of the satellite, and to obtain good approximations for it. 
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Because it seems to be impossible to get an explicit solution for the eleva- 
tion angle error ♦ a computer program (which will be documented in a com- 
panion report) was designed to obtain accurate numerical solutions. This 
program can also evaluate any approximate formula for , and may be used to 
test the accuracy of such formulas. The approximate formula must be coded in 
a standard form and added to the main deck. The refractlvlty profile (l.e., the 
variation of the index of refraction with altitude) is arbitrary; it also is coded 
in a standard form and added to the main deck. The program as used to test 
the accuracy of the various approximations obtaiiied for , and some of them 
were found to be quite accurate. 

The work presented in this report was inspired by the writer’s desire to 
have a mathematically rigorous, systematic and detailed treatment of the prob- 
lem of the refraction of a ray through a spherically symmetric atmosphere. It 
is felt that the treatment of the problem given below satisfies those require- 
ments reasonably well. Many of the results appear elsewhere in the literature, 
and some are well-known. Some of the results, however, are believed to be new; 
e.g., the simple upper bound for the elevation angle error given by the inequality 
(30), the very genera] approximate formula for the same quantity given by Equa- 
tion (61), and the inequality (23), which gives a necessary and sufficient condition 
for the existence of a ray path joining a satellite to a ground station, in terms of 
the satellite’s position and the refractivity profile. Even though some of the re- 
sults obtained are well-known, the writer feels that the very completeness of the 
mathematical treatment given in this report is valuable in itself, because it 
brings to light questions (such as the existence of ray paths and their general 
behavior) which are usually not consldertvi. 


BACKGROUND IklATERIAL 

The material in this section is not necessary for the mathematical discussion 
which follows. It is inserced here to provide some background information and 
to motivate the various assumptions made in the mathematical treatment. 


Tropospheric Index of Refraction 


In most regions of the electromagnetic spectrum the troposphere is es- 
sentially non-dispersive; i.e., the index of refraction does not change appreciably 
with frequency. The index of refraction n of a medium is defined by n = c/v, 
where c is the velocity of light in vacuum and v is the velocity of light in the 
medium. Because the value of n is nearly unity throughout the atmosphere (in- 
cluding the ionosphere) , it is customary to use a related quantity called the 
’’refractivity," defined by N = (n - 1) x 1C®. The following empirical formula 
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for the radio refractivity of the troposphere is in widespread use by radio 
scientists (see Ref. 2, p. 7): 




4810 e 


( 1 ) 


where T is the temperature ii» degrees Kelvin, P is the total atmospheric pres- 
sure in millibars and e is the partial pressure of water vapor in millibars. This 
formula is considered to be accurate to 0.3% for radio frequencies up to 30,000 
megahertz, for normally encountered ranges of pressure, temperature and 
humidity. The fact that the formula does not involve the frequency shows how 
nearly frequency-independent the refractivity is in this region of the spectrum. 

The ratio e/P is about 0.01 at ground level for 60% relative humidity, and 
decreases rapidly with altitude (Ref. 2, p. 11). It follows that the second term 
in Equation (1) is at least an order of magnitude smaller than the first term. 
Since the pressure P (which is the dominant term) normally decreases expo- 
nentially with height, and since the variation of T with height is not great, it 
would seem that an exponential decrease of N with height might be a good ap- 
proximation to the actual refractivity profile. In Ref, 2 it is shown that an ex- 
ponential radio refractivity profile is in fact a quite accurate model for the 
troposphere. 

For optical frequencies the index of refraction of air is also very weakly 
frequency-dependent. The refractive bendii^ of violet light differs from that of 
red light by only about one percent. The dependence of N on pressure and tem- 
perature, for visible light, is given in Ref. 3 (p. 62) by the formula 


p 

N = 82- 
T 


( 2 ) 


where N, P and T are as defined above (the units used in Ref. 3 have been 
changed to the units used here) . This formula is almost the same as formula 
(1) , with the second term deleted. It thus appears that the tropospheric refrac- 
tivity profile for visible light should also be closely approximated by an ex- 
ponential model. 
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Ionospheric Index of Refraction 


The refractive behavior of the ionosphere is very different from that of the 
troposphere. The density of air above 40 kilometers is too small to cause any 
appreciable refractive effects; the refractive behavior of the ionosphere is due 
primarily to the rarified electron gas arising from ionization by solar ultra- 
violet and X-radiation. (For a recent review of what is known about the forma- 
tion and composition of the ionosphere, see Ref. 4.) A discussion of the propa- 
gation of electromagnetic waves in such a medium is given in Ref. 5. The follow- 
ing formula is derived there (p. 329) for the index of refraction n of an electron 
gas in the presence of a static magnetic field (the earth's field): 


Me 2 
me„ 


1 - 


i 


e/Xo Hq o) 


m 


(3) 


where M is the number of electrons per unit volume, e is the electronic charge, 
m is the electronic mass, is the permittivity of free space, Mq is the per- 
meability of free space, is the intensity of the component of the magnetic field 
in the direction of propagation of the wave and a> is the circular frequency of the 
wave. All quantities are in the rationalized m.k.s. system. This formula shows 
that the ionosphere has some remarkable refractive properties. Perhaps the 
most striking thing about the formula is the presence of the ± sign; this means 
that there are two modes of propagation with two distinct velocities . In fact it 
can be shown that a linearly polarized wave entering the medium is resolved 
into right- and left-circular ly polarized components; thus a rarified electron 
atmosphere with a static magnetic field imposed on it acts like an anisotropic, 
doubly-refracting crystal. Further, it is apparent from Equation (3) that for a 
certain frequency the index of refraction of one of the waves becomes infinite, 
and that for another frequency the index of refraction of the other wave becomes 
zero. This behavior accounts for many of the remarkable properties of the 
Kennelly-Heaviside layers at radio frequencies. 

The effect of the earth’s magnetic field becomes less important for higher 
frequencies, as can be seen from Equation (3). In fact, if we use the values 
e = 1.602 X 10~ coulomb, e/m = 1.759 x lo“ coulombs/kilogram, = 4rr x 
10'^ henrj'/meter and = 8.854 x 10~^2 farad/meter, and if we take the maxi- 
mum value of to be about 0.4 gauss, then we see that the ratio of the two 
terms involving o) in Equation (3) is 

1.407 
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where f — (jO /2 Tf is the frequency in cycles per second. Hence for f > 100 mega- 
hertz, which is about the lowest frequency likely to be used for satellite track- 
ing, the above ratio is less than 1.4 x 10”'*. We may therefore neglect the effect 
of the earth's magnetic field and represent the ionospheric index of refraction 
by the simpler formula 


,2 .. 


1 - 


Me^ 


m^o 


(4) 


This simplified formula still has some interesting properties. For the critical 
frequency w = = (Me^/meg)* ^ the index of refraction becomes zero. Physi- 

cally this means that a wave entering the medium is totally reflected. More 
generally, if M is increasing monotonically along the direction of propagation, a 
wave of frequency cj will be totally reflected when it reaches the region where 
M = me^ oj^/e^. This phenomenon is used to measure the electron density pro- 
file of the ionosphere by beaming radio pulses vertically upward and measuring 
the time between transmission and reception of the pulses. As the transmitter 
frequency is increased, higher regions of the ionosphere are probed. Only that 
part of the ionosphere below the F 2 majdmum (which lies between 200 and 400 
kilometers) can be probed by this method, however, because the electron density 
above that point decreases again. To measure the top part of the ionosphere it 
is necessary to use satellites to beam radio pulses downward. 

If we use the values given above for e, m and , and if we assume that the 
largest value of M in the ionosphere is about 10 * ^ electrons per cubic meter (a 
realistic assumption, corresponding to the F^ maximum), we find that =<^o / 
2 t 7 = 8.98 megahertz. This is far below our lower bound of 100 megahertz; we 
therefore do not need to concern ourselves with the possibility of total reflection. 

There is an interesting property of the above formulas for n which we have 
not yet mentioned: the value of n is less than unity. This means that a wave of 
any frequency greater than f,. is propagated in a rarified electron gas with a 
velocity greater than the \ ciocity of light in free space. This of course does not 
contradict relativity theory. In the derivation of Equation (3) it is implicitly as- 
sumed that a sinusoidal plane wave of frequency f, of infinite extent, is beii^ 
propagated through a rarified electron gas which is also of infinite extent. But 
every physical wave is limited in space and time, and has a frequency spectrum 
of finite (i.e., non-zero) width. Such a waveform can be represented mathemati- 
calh y a Fourier integral, which may be regarded intuitively as the sum of an 
(uncountably) infinite number of sinusoids of infinitesimal amplitude, each sinu- 
soid being of infinite extent. Each of these hypothetical sinusoidal waves is 
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propagated with a velocity greater than c, the velocity of light in free space; but 
their sum, the actual physical waveform, is propagated with a ’Velocity” less 
than c. The reason for using quotation marks is the impossibility of precisely 
and unambiguously defining the velocity of the physical wave. Each of the com- 
ponent sinusoids is propagated with a different velocity, because the index of 
refraction varies with frequency (i.e., the medium is ’’dispersive”). Therefore 
their phases are constantly varying relative to each other, and the shape of their 
sum is constantly changing. The physical waveform is thus continuously deformed 
as it moves through the medium, and it is not possible to uniquely define its 
velocity. However, it is possible to define velocities which are useful under cer- 
tain circumstances. The one most commonly used is the group velocity; it is 
meaningful only for waveforms whose frequency spectrum is confined to a very 
narrow band. If this band lies in a region where the dispersion of the medium is 
normal and moderate (e.g,, away from absorption bands), the group velocity is 
approximately the velocity of energy propagation. In Ref. 5 the following expres- 
sion is derived for the group velocity: 


= -L 

''group ~ 

dco 


where k is the wave number (k 277 /A., where is the wavelength in the medium), 
and the derivative is evaluated at , the center of the narrow frequency band. 

If we use the relation k = a> /Vphase (oj), where (cy) is the velocity (in the 

medium) of the component sinusoid of circular frequency o), we obtain a relation 
between the group index of refraction n (= c/v ) and the phase index 

* group ' ' group ' ^ 

n . (= c/v ^ ): 

phase ' phase ' 


n , " n , 

group phase 


1 + 


d n 


phase 


n 


phase 


do) 


(5) 


The phase velocity is so called because it is the rate at which surfaces of 
constant phase, of a pure sinusoidal plane wave of infinite extent, are propagated 
through the medium. If the medium is non-disp)ersive (as the troposphere is, for 
example), the derivative is zero and the group velocity is the same as the phase 
velocity (which is then independent of frequency). 

If we apply Equation (5) to the expression for given by Equation (4), 

we obtain the following expression for the group index of refraction of the iono- 
sphere (neglecting the earth's magnetic field); 
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n = 1 - 

group 1 


M 


-1 2 


( 6 ) 


m 


It is apparent from Equation (6) that v < c, as it should be, if the group 
velocity represents the rate at which energy is propagated through the medium. 

It is also worth pointing out that v , • v . = . 

^ ^ group phase 

The problem of the propagation of electromagnetic signals through a dis- 
persive medium is a very complex one. A very illuminating discussion of this 
problem is given in Ref. 5; it is based on an investigation conducted by Sommer- 
feld and Brillouin in 1914 (Ref. 6) The results of this investigation show that the 
time t of arrival of a signal at a given point in a dispersive (homogeneous and 
isotropic) medium can never be less than d/c, where d is the distance of the 
given point from the signal source, which is assumed to begin emitting the sig- 
nal at time t = 0. 

It is important to note that the phase index must be used to compute the 
bending of a radio wave through the ionosphere, but the group index must be used 
to compute the travel time of the signal. To see why this is so, consider a signal 
which has a narrow frequency spectrum, in the sense that it has a small ratio of 
spectral width to carrier frequency, (A modulated radio signal will ordinarily 
have a narrow spectrum in this sense, and even a radar pulse consists of so many 
cycles that its spectrum is fairly narrow.) Each (hypothetical) component wave 
of frequency f siaffers refraction, the degree of which depends on its index of re- 
fraction n(f), which we have called the phase index. Because of the narrowness 
of the spectrum, the variation (throughout the spectral width) in the amount of 
refraction is small. In other words, the angular dispersion of the signal in space 
is small. Thus the curved path followed by the physical wave is essentially the 
same as that followed by any one of its hypothetical components, for the calcula- 
tion of which the phase index must naturally be used. 

On the other hand, if the travel time is to be computed the velocities of the 
hypothetical component waves are not of interest. It is the velocity of the re- 
sultant physical wave that matters; and this, as we have seen, is given by the 
group index of refraction. 

Good approximations for the phase and group refractive indices may be 
derived from equations (4) and (6) by using the binomial approximation (1 - a)*^ 

1 ^ a/2, with a = Me^/megO;^. If we use the values for the constants given 
previously, and if we again take 10^^ electrons per cubic meter as the maxi- 
mum value of M and 10® cycles per second as the minimum value of f, we see 
that the maximum value of a is about 0.00806. The binomial approximation is 
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thus quite accurate. If we use it, and if we insert the values of the constants, 
equations (4) and (6) become 


and 


n 


phase 


= 1 


40.3 iL 
f2 


n = 1 

group 


+ 40. 3 ^ , 


(7) 

( 8 ) 


where M is the electron density in electrons per cubic meter and f is the carrier 
frequency in cycles per second. In terms of the refractivity N = (n - 1) x lO® 
these equations become 


and 


N 


phase 


- 40.3 X 10® ii 
f2 


N 

group 


40.3 X 10 ^ — . 

f2 


(9) 


( 10 ) 


From these formulas it is clear that both the phase and group refractivity pro- 
files in the ionosphere are proportional to the electron density profile, for a 
given frequency, and that they decrease rapidly with increasing frequency. For 
radio frequencies of a few hundred megahertz the ionospheric contribution to 
refraction is important. For frequencies of thousands of megahertz it becomes 
very small, and for optical frequencies it is infinitesimal compared to the ef- 
fects of the troposphere. 


The Approximation of Geometric Optics 


In many published papers dealing with ionospheric refraction the analysis 
is based on the laws of geometric optics, with no comment as to the validity of 
these laws in the ionosphere. It is not difficult to derive the laws of geometric 
optics from Maxwell’s equations, for the case of dielectric media; such a deriva- 
tion is given in chapter 3 of Ref. 1. The derivation is more complicated for a 
plasma such as the ionosphere. It turns out that the laws of geometric optics do 
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hold for the ionosphere, provided the spatial variation of electron density is suf- 
ficiently gradual; just as for dielectrics, these laws will hold provided the prop- 
erties of the medium do not change much over a distance of one wavelength of 
the radiation. Fortunately this condition is usually satisfied for radio waves in 
the ionosphere. Very detailed analyses of the applicability of geometric optics 
to ionospheric propagation are given in References 7 and 8. 

The main result from geometric optics that will be used in the ensuing 
analysis is the spherical form of Snell’s law, sometimes called Bougiier's Law. 
It applies to a spherically symmetric medium, and states that 


n R sin 6 = constant , 


where n is the index of refraction at the distance R from the center of symmetry 
(in our case, the center of the earth), and 6 is the angle between the ray path 
and the radius vector (from the center of symmetry) at distance R from the 
center. This law is the spherical analogue of the more commonly used planar 
form of Snell’s law. It should be mentioned that the index of refraction n is the 
phase index, for the frequency of the radiation being used. 


FORMULATION AND SOLUTION OF THE PROBLEM 
OF FINDING RAY PATHS 

Formulation of Problem 


Brief Statement of Problem - We want to solve the following problem: 

Given a point source of radio or radar waves in or above the earth’s atmosphere, 
and given a station on the ground, does there exist a ray path joining the source 
to the ground station? If so, how can we compute the path? In particular, how 
can we find the apparent elevation of the source as seen from the ground? The 
ray path will of course deviate from the slant path, and the measured elevation 
angle will differ from the true elevation angle, because of atmospheric refrac- 
tion of the radio waves. 

Simplifying Assumptions — To make the mathematical analysis tractable 
we make the following assumptions: 

• The atmosphere is spherically symmetric; i.e., the index of refraction 
depends only on height (above a spherical earth). 

• The effects of the earth’s magnetic field are negligible. 
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• The propagation of radio waves through the atmosphere can be adequately 
described by the laws of geometric optics (i.e., by ray-tracing methods). 

The second and third assumptions will generally be valid for radio frequencies 
above 100 megahertz. The first assumption is only a rough approximation, since 
the ionospheric refractivity profile varies considerably with latitude. 

Boug uer's Law as the Characterizing Property of Ray Paths - With the 
above assumptions we can briefly state our problem as follows: Find a curve 
(the ray path) joining the source to the ground station, which has the property 
that it satisfies the spherical form of SnelPs law (Bouguer’s law). 


n R X s = constant , 


( 11 ) 


at every point along the curve. (See Chapter III of Ref. 1 for a derivation of this 
relation.) In particular, find the arrival angle; i.e., the angle between the ray 
path and the horizontal at the ground station. Figure 1 illustrates the geometry. 
In equation (11) and Figure 1, R is the radius vector from the center of the earth 
to any point on the ray path, R is the magnitude of R, n is the phase index of re- 
fraction (for the radio frequency being used) at distance R from the center of the 
earth and s is a imit vector tangent to the ray path, directed outward along the 
path (i.e., away from the ground station). The arrival angle is denoted by , 
and the true elevation angle of the satellite by a^. The angle 4 ^ is the central 
angle, as shown in Figure 1. The other quantities are defined below. It is shown 
in References 7 and 8 that the derivative dn/dR must be continuous, and in fact 
slowly varying, in order for geometric optics to be applicable. We will there- 
fore assume that n is a continuously differentiable function of R, from R^ (the 
radius of the earth) to infinity, and we will look for ray paths which are continu- 
ously differentiable curves. 

Planarity of the Ray Path - It follows from Equation (11) that any ray path 
(i.e., any curve which satisfies Equation (11)), lies entirely in a plane containing 
the earth’s center. To see this, consider any point p along the curve, and let 
Pp be the plane determined by the vectors R and s at the point p. N(w let q be 
any other point along the curve, and let P^ be the p&,ne determined by R^^ and 
^ . It is clear from Equation (11) that the normals to P^ and P^ are parallel; 
since both planes must pass through the center of the earth, it follows that they 
must be identical. Since p and q are arbitrary points, it is clear that the entire 
curve must lie in one plane. 
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Figure 1-The Ray Path 


Boiler’s Law in Scalar Form - If there exists a ray path which joins the 
satellite to the ground station, then of course it lies in the plane determined by 
those two points and the earth's center. The ecpiation of such a ray path can 
therefore be written as a single scalar equation. Probably the most nabiral way 
of deriving a scalar equation from Equation (11) is to take the magnitude of each 
side; if we take account of algebraic signs, the result can be written in the form 


n R sin 6 = c, 
1 


( 12 ) 


11 


wher^i^ is the angle between R and s , taken positive in the clockwise sense 
from R toward s , with ~v < 6 and c is a constant. Equation (12) applies to 
any ray path, whether or not it passes through the ground station or the satellite. 
If we explicitly require that the ray path terminate at the ground station by evalu- 
ating the constant c at that point, Equation (12) becomes 


n R sin 6 = n R sin 0 
g c g 


= n R, cos a 

g e 


a ’ 


(13) 


where n , (the radius of the earth) and 6^ are the values of n, R and 0 , re- 
spectively, at the ground station and = v/2 - Q . Note that Equation (13) is 
merely a condition that a ray path must satisfy; it does not explicitly give the 
ray path in any coordinate system. In fact, Equation (13) is really a differential 
equation for the ray path, because the angle 6 is clearly related to the slope of 
the curve. 

The Epsilon Function - It is convenient to introduce a quantity e, related to 
the index of refraction n, and defined by the relation 



n 


(14) 


Since n differs only slightly from unity, it follows that e is a small quantity; 
further, e (R^ ) - 0, since n(Rg) = n^. It is also worth noting that the condition 
0 < n (R) < oo (for all R > R^ ) implies that _ l‘^e(R)< °° (for all R > R«); this 
condition will certainly hold in the earth's atmosphere, including the ionosphere, 
for frequencies above 100 megahertz. In terms of e , Equation (13) becomes 


sin ^ r (1 + e) 



cos a . 

- a 


(15) 


Equation (15) is the relation which must be satisfied by a ray path emanating 
from the ground station at the angle . 

Condition for the Existence of a Complete Set of Ray Paths - It is important 
to note that there may not be any curve which satisfies Equation (15), for a given 
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value of o.^ . For such a curve to exist, it is obviously necessary that the right- 
hand side of Equation (15) be less than or equal to unity for every value of R in 
some interval starting at . If we want to ensure the existence of a ’’complete 
set” of ray paths, l.e., a ray path of infinite extent (emanating from the ground 
station) for every value of in the interval 0 < 5 tt/2 , it is necessary to as- 

sume that the condition 



is satisfied for every value of R in the Interval R^ < R < oo. If we solve this in- 
equality for £ , it becomes 


e (R) < 


R -R. 


for all R ^ R^ . It turns out that this condition is not sufficient for the existence 
of a complete set of ray paths, but it will be shown that the slightly stronger 
condition 


R _ R 

- 1 <e (R) < (1-a) ♦-_! 

K 


(16) 


for all R ^ R^ , where 0 < a < 1 and equality holds only for R = R^, is sufficient, 
(The fact that e > - 1 is obvious from the definition of £ .) 

The Class £ of Epsilon Functions - For convenience we will denote by 8 
the set of e -functions which have the following properties: 

• e (R) is continuously differentiable for all R ^ R^ ; 

• € (R) satisfies condition (16) for some value of cr in the interval (0,1). 

It is worth noting that it is not clear whether the condition (16) is necessary for 
the existence of a complete set of ray paths, and it seems possible that it night 
ever i -2 unduly restrictive. It turns out that it is not too restrictive, at least in 
the sense that any physically realistic profile n(R) will probably satisfy the 
conc'ition (see the appendix for an example). 
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Derivation of an Explicit Equation for Ray Paths 


We will now show that, for any e-function In the class 8 , equation (15) can 
be written as a differential equation In the variables R and >//, for any a^e[0, tt/2), 
and that th^ s equation can be integrated to yield an explicit equation for the ray 
paths (with as a parameter). 

Some General Properties of Ray Paths — We will begin the proof by showing 
that, for any e [0, 77^2), the ray path (assuming it exists) will have the property 
that R is a single-valued and strictly monotonically increasing function ot 4>. It 
is clear that such a path will never dip below the earth's surface and that it will 
never "double back" on itself across any given radius vector. 

The case = tt/ 2 must be treated separately. It is clear that In this case 
the ray path will just be a vertical straight line, whether or not condition (16) is 
satisfied, because such a line satisfies Equation (15). No other path can be a ray 
nath in this case, since such a path (which necessarily would be initially vertical 
and would later deviate from the vertical) would have to intersect some non- 
vertical line from the center of the earth; at the point of intersection, sin d would 
not be zero. (It is also clear from the symmetry of the situation that the ray 
path must be a vertical straight line.) 

We now return to the general case where < tt/2 . We first show that each 
ray path has the property that R is a single-valued function of i whether or not 
condition (16) is satisfied. From Equation (15) it is clear that, for any e [0, 77/2), 
sin 5^ > 0 at every point along the associated ray path, assuming one exists. 

Since -77 < + 77 (by definition), it follows that 0 < ^ < 77 at every point along 

the path. It is therefore clear that the ray path cannot double back across any 
given radius vector; in other words, for each value of 0 there is at most one 
value of R. 

We next show that, for each e [0, 77/ 2), the associated ray path (if It 
exists) has the property that R is a strictly monotonically increasing function 
oi4> t if condition (16) is satisfied. If we use the inequality (16), it follows from 
Equation (15) that 


s in t? 5 



' cos a f cos a 

m a 


for all R > R^, with equality holding only for R = R^. It follows that sin 8 < 1 
at all points of every ray path, except that sin <9 =1 at the single point R = R^ 
for the ray path for which = 0. We recall that we are lookii^ for ray patl^ 
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which are continuously differentiable curves, and are therefore continuous; for 
such a ray path it follows from the above condition on sin 6 that 6 < v/2 at every 
point, except of course that 6= v/2 at the ground station for the ray path for 
which = 0. It follows that for each such ra’’ ^>ath R increases strictly mono- 
tonically with . 

Transformation of Bouguer*s Law to a Differential Equation in Polar Co- 
ordinates — We are now ready to write Equation (15) as a differential equation. 
This will be accomplished by deriving a relation between the angle 0 and the 
derivative d>p/dR; the desired relation can be easily found with the help of 
Figure 2. 



Figure 2-Diagram for the Calculation of dR/6\p 




From the figure we see that 



and 


tan 6 = 

dR 


The second equation is invalid only for the point R = along the ray path (if it 
exists) for which = 0. With this one exception, it follows that 


dR _ R 
d\p tan 6 


The presence of tan 8 in the denominator causes no difficulty; as we have shown 
above, 0 < ^ < v/2 along every ray path, and therefore 0 < tan 0 < co (except that 
tan 6 is infinite at the one point mentioned above). It follows that 0 < dR/d^ < oo , 
except that dR/d0 = 0 at the one exceptional point. Because of the strict mono- 
tonicity of R(0) we know that the inverse function 0(R) exists; from the above 
condition on the derivative of R(^) we have 0 < di/^/dR < qd , except that d</'/dR = oo 
at the exceptional point. We may therefore write 


dip _ 1 _ tan 6 
dl'' ^ R 
d0 


From Equation (15) we see that 


(17) 


tan 8 = 


sin 6 


+ /l - sin^ 8 


(1 + e ) * — • cos a 
R ° 

' T" 

1 - (1 + OM — ) cos^a^ 


(18) 
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The plus sign is included to explicitly show that cos 6> 0, since 6 is known to 
be in the first quadrant; this equation is therefore completely equivalent to 
Equation (15). From Equations (17) and (18) it follows that 


d0 

dH 


( 1 . 0 - 



• cos a 

a 


-(1 + e )- 



cos^ a 

a 


(19) 


This is a differential equation for the ray path, for any given angle a^. The 
quantity under the square-root sign will always be greater than zero, except at 
the one exceptional point, where it is zero. (It should be noted that this may fail 
to be true if condition (16) does not hold.) 

The Explicit Equation of a Ray Path — We have shown that, for any e e 8 , the 
definii^ equation for a ray path may be written in the form of the differential 
equation (19). It must now be asked whether this equation has solutions; if it 
does, then of course those solutions are the ray paths we are looking for, and 
they have all the properties described above (such as monotonicity). To find 
out whether Equation (19) has solutions for all we integrate it formally to 
obtain the relation 



[ 1 + e (s)] -L cos a 
2 “ 


/R 

1 - [1 + e(s)]^ ( — j cos^ 


d s 


(20) 


which is the equation of a ray path in the polar coordinates \p and R , provided the 
integral exists. For 0 the quantity in the denominator is greater than zero 
for all s > , and the integral therefore exists. The case = 0 requires a 

special proof of the existence of the integral, since the integrand becomes in- 
finite at R = R^ . The proof makes use of the inequality (16) to obtain a majorant 
function for the integrand; the integral obtained by replacing the original inte- 
grand by this majorant function is then shown to exist. If the condition (16) is 
assumed to hold, then a simple calculation shows that 
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1.(1. o' 



for all s > , with equality holding only for s = R^ . It also follows from (16) 

that 


(1 + e) 1 -C7 +(T 



1 

s 


< 


^ 1.0 + 1 


Re / ^ • '/i' 


2 

• /s 

e 


for all s > R^ . With these inequalities we may form the majorant function M(s): 

' ^ 2 


(l.e) 


yi-(1.0' (^) 


R \ 2 /T • /a (S-RO 


= M(s), 


for all s > R^ . As s -* R^ , both functions become infinite. It is easy to show that 
the integral exists with this majorant function as the integrand; the substitution 
X = s - R^ yields 



M(s) ds 




- 1 . 
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which is finite for all R > . It follows that the integral exists even when 

= 0. (It is not difficult to show that if we were to allow cr to be zero the inte- 
gral would diverge for = 0.) We conclude that, for any € e 8 , a complete set 
of ray paths exists; they are given explicitly by Equation (20). 


The Existence of a Ray Path Connecting a Given 
Satellite Position to the Ground Station 


We will now investigate more closely the character of this bundle of ray 
paths. First, we note that no two ray paths can intersect. To see this, it is 
sufficient to note that the integral in Equation (20) is a strictly monotonically 
increasing function of cos ; it therefore cannot take the same value (for a 
given value of R) for two different values of . Because no two paths intersect 
it is clear that the ray path (if it exists) joining a given point to the ground sta- 
tion is unique. 

The second property worth noting about the bundle of ray paths is that the 
paths change continuously with . A rigorous statement of this property- would 
require that the set of ray paths be made into a normed space, perhaps with the 
"sup norm;" but the property can be seen intuitively by noting the fact that the 
integral in Equation (20) is a continuous function of , for any given value of R. 

The third and most important property of the set of ray paths concerns the 
region of space which is "covered" by the paths; i.e., the set of points which lie 
on some ray path. To investigate this question will be our next task. 

General Condition for the Existence of a Ray Path - Suppose we are given a 
point (a^, R^), where is the (true) elevation angle of a satellite (taken positive 
upwards), and R, is its distance from the earth’s center. We can determine 
whether this point Ues on one of our ray paths by setting 4 > ~ 4 >^ (f^® value of <// 
corresponding to and R ,) and R = R, in Equation (20); if there exists a value 
of in the interval [ 0,7t/2) which satisfies the resulting relation, then the point 
lies on the ray path corresponding to that value of . 

To write the relation explicitly we must express in terms of and R^ ; 
this is easily done by applying the law of sines to the triangle formed by the 
earth’s center, the ground station and the satellite, as shown in Figure 3. The 
law of sines states that 


cos 


s in 7 
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SATELLITE 



CENTER OF EARTH 


Figure 3-Diogram for the Calculation of 


or 


Re 

Sin y - — cos . 
R ‘ 


This expression for sin y is easily seen to be valid for and e (-7t/2, 

71 ^ 2 ) ; i.e., for any point outside the earth and to the right of the ground station. 
It is clear that sin 7 < 1 everywhere in this region; since 7 is clearly a continu- 
ous function of and R^ in this region, and since it is obvious that 7 < tt/2 
for some points, it follows that y < 2 at every point of the region. We may 

therefore write 


7=sin"^ — cos a. 

\ R. ‘ 
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where the value of the arc sin is taken in the first quadrant. It follows that 


= 7T - 

~s 



- y 


or 



( 21 ) 


where the value of the arc cos is taken in the first quadrant; this expression for 
is valid for R >R and a e (-77/2, 77/2) . 

s set' ' 

We can now state that the point (a^ , R ) will lie on some ray path if and only 
if the relation 



R, 

(1 + e) — cos dR 


1 * (1 I cos2 a 


1 2 



(22) 


is satisfied for some value of a.^ £ [0, 77 ^2). Since the integral decreases mono- 
tonically as increases, we see that Equation (22) will have a (unique) solution 
if and only if 


r 


(l + 6)_dR 
R2 


/ R ^ 2 


1/2 


^ cos*^ 




■ cos . 


(23) 


In other words, given any e -function in the class S and any point (a^ , R,) ^ 

77 /2) X (R^ , 0 D), there will exist a ray path connecting that point to the ground 
station if and only if condition (23) is satisfied. 

We have now answered the question posed earlier concerning the region of 
space covered by the complete set of ray paths (which is guaranteed to exist for 
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any e -function in 8 ) ; we have shown that, for any given e -function in 8 , the points 
which lie on some ray path are just those points which satisfy condition (23). As 
it stands, this condition is not very transparent; to get explicit results we must 
consider some special cases. 

The Existence of Ray Paths in Some Special Cases — One special case is 
suggested by the practical requirement that the satellite be on or above the 
horizon; we might therefore be interested in values of greater than zero. 

Since the expression on the right-hand side of (23), considered as a function of 
, decreases (strictly) monotonically for e L-tt/Z 2], we see that (23) will 
hold for all e [0, tt/ 2 ] if and only if 



(1 +e) _ldR 
R2 


1 - (1.C)2 



1/2 


« COS“ ^ 



If we note that 



Re 

(1 + €) — dR 
R2 


1 - ( 1 + 6)2 



2 



> 



= cos“* 



(24) 


for any e -furecion in 8 which has the property that e (R) > 0 for all R, we see 
that condition (23) does in fact hold for all a^e [0, n/ 2] and for all R, , for any 
such € -fimction. In other words, if e (R) is everywhere non-negative, then the 
entire region in and above the horizontal plane and to the right of the ground 
station is covered by the complete set of ray paths. Of course the condition e (R) 
^ 0 is merely sufficient, not necessary, for this region to be covered. Also, for 
a given non -negative e -function the total region covered will be even larger, 
unless e (R) = 0. 
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To see the truth of this last statement, suppose that e (R) ^ 0 for all R, but 
e (R) 0; then (24) becomes a strict inequality. From the fact that the right- 

hand side of (23) decreases monotonically with for -tt/2 ^ < 77/2, it then 

follows that, for each value of R^ , there exists / 3 e [-77/ 2 , 0) such that 



for all > yS , for the given value of R , . Comparison of this relation with (23) 
shows that for any non-negative £ -function which is not identically zero, and for 
any value of R, , there is an angle /3 < 0 such that every point (a^ , R^) for which 
> /9 is connected to the groimd station by a ray path. If, in addition, it happens 
that 


/3„ = sup {^(R.):R. >R.><0, j 

I 

then the entire region I 

= An - ^ 

is covered by ray paths. 

The condition € (R) > 0 will actually be satisfied by a realistic refractivity 
profile, provided that unusual conditions, such as temperature inversions near j 

the ground, do not occur. To see this, one might note that the tropospheric re- | 

fractivity is commonly assumed to decrease exponentially with altitude; thus 
n < n , and therefore €=n^/n-l>0, for all R > R ^ . In the ionosphere n < 1 
and therefore (since n^ > 1) we again have e > 0. 

As another special case, suppose we ask what happens if € (R) 1 0 for all R, 
but e (R) * 0. One thing at least is easy to see in this case: no point on the 
horizontal plane lies on a ray path. To see this we just note that for such a point 
to lie on a re^ path the condition 


i 
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R 

(1 +e) -IdR 
R2 


1 -( 1 + 6)2 



2 

_ 


1/2 


^ COS 


1 



must be satisfied; but for 6 (R) ^ 0 and not identically zero this is impossible, 
because for such an e -function 



Since no ray path passes through or even touches the horizontal plane, it follows 
that all ray paths lie completely above that plane for any non-positive (and non- 
zero) € -function. 

More detailed results on the existence of ray paths could be derived, but for 
our purposes we do not need to go any further with this analysis. 


Recapitulation 


Summary of Results - Before continuing with our analysis we will recapitu- 
late what we have found so far. We showed that any ray path jc.'ning a satellite 
to a ground station throi^h a spherically symmetric atmosphere must lie entirely 
in the plane determined by the satellite, the ground station and the center of the 
earth. Since there may not be ftny ray path joining a given satellite position to 
the ground station, it seemed most logical to focus our attention on the bundle of 
ray paths emanating from the ground station, rather than on the satellite position; 
we may limit this bundle to a plane, because of symmetry, and may {uirametrize 
it by the variable , the elevation angle of a ray path at the ground station. We 
stated a condition on the refractivity profile (condition (16)) which we claimed is 
sufficient to ensure the existence of a complete set of ray paths; i.e., a ray path 
of infinite extent for every value of in the interval [0,7r/2l. (This condition 
is actually satisfied by realistic refractivity profiles; see the appendix for an 
example.) On the assumption that condition (16) holds, we then actually proved 
that a complete set of ray paths exists; we derived an explicit equation for them 
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(Equation (20)), with as a parameter, and we showed that they have certain 
properties. Among these properties are a single -valued dependence of R on 0 , 
with R increasing strictly monotonically with </' (for < tt/2). We also showed 
that no two ray paths intersect, which is equivalent to saying that the ray path 
joining a satellite position to the ground station is unique, if it exists at all. We 
next Investigated the question of the region covered by the bundle of ray paths. 

We derived a condition (condition (23)) which is both necessary and sufficient 
for the existence of a ray path joining a given satellite position to the grountl 
station, for a given refractivity profile. We then investigated some special cases; 
the most important of these, and the one v© will be exclusively concerned with 
in the rest of this report, is the case where the e -function is non-negative for 
all R. In this case we showed that any point above the horizontal plane is joined 
to the ground station by a ray path. This fact, together with the fact that real- 
istic refractivity profiles actually have the property that e (R) > 0 for all R, 
justifies concentrating our attention on this type of profile. 

Procedure for Calculating a Ray Path - We can now outline the procedure 
for finding a ray path. Suppose we are given an e -function (in the class P ) and a 
satellite position; to determine whether a ray path exists, we test to see whether 
these satisfy the inequality (23). If the inequality is not satisfied, there is no ray 
path; if it is satisfied, there is a unique ray path joining the satellite position to 
the ground station. If e (R) > 0 for all R and the satellite position is above the 
horizon, we know that the inequality will be satisfied. To find the ray path, if it 
exists, we must first solve equation (22) for ; this value of can then be put 
into equation (20) to obtain the equation of the ray path in the polar coordinates 
R and 4i . The only step in this procedure which offers any difficulty is the solu- 
tion of Equation (22); unfortunately this equation cannot be explicitly solved for 
, and we must resort to a numerical solution or to approximate formulas for 
in terms of and R, (which describe the satellite position). The rest of 
this report will 1^ mainly devoted to finding such approximate solutions and 
numerical solution methods; we will also try to deduce, directly from Equation 
(22), as much information as possible concerning the properties of the fimction 

“a(^. R,)- 

INVESTIGATION OF THE ELEVATION ANGLE ERROR AS A FUNCTION OF 
SATELLITE POSITION, FOR A SPECIAL CLASS OF REFRACTIVITY PROFILES 

The Class 8 ^ of Epsilon Functions 


As was mentioned above, we will limit our attention in the rest of this re- 
port to e -functions which are non-negative, and to satellite positions on or above 
the horizon (a^ > 0). 
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For convenierce we will denote the class of non-negative, bounded e -functions 
by P* ; l.e., 

c* = {e I € e P, e (R) > 0 for all R - and € (R) is bounded on [R^, <»)} • 

(Note that P •) I*' other words, an e -function is in the class £ if and only if 
it has the following properties: 

• 6 (R) is continuously differentiable for all R ^ ; 

• e(R)^ 0 for all R> R^ ; 

• There exist a number cr e (0,1) and a positive number B such that 

f (R) < Min |(1 -cr) • 

for all R > R^ . 

We will also denote by D the region of the - R^ plane in which we are 
interested: 

D = {(a^ . RJ 1 0 < < tt/2 and R, > R.) 

= [0, 77/2) X (R,.®). 

A point (a^, R^ ) In the - R, plane will sometimes be denoted by p, for brevity. 

The Equation for the Elevation Angle Error 

For > 0 it is possible to write Equation (22) in r. more transparent form, 
by expressing the right-hand side as an integral similar iti to the integral 
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on the left-hand side. We only need to note that, for e (R) - 0, the ray path will 
be a straight line and ; we can therefore use Equation (21) and Equati^*- 

(20), with e (s) =0, a, = aj and R = , to write 



The correctness of this expression can be easily checked by direct evaluation of 
the integral, using the substitution (R^ /R) cos = cos x. Thus, for > 0, 
Equation (22) takes the form 



(1 + €) — cos (a^ + 8g) dR 


il - (1 e)2 ( ^ j cos2(a^+Sg) 



R„ 

— cos a. dE 
R2 



= cos 



- a 


t ’ 


(26) 


where we have written +8 . The quantity 8 is the "elevation angle 

error;" i.e., it is the difference between the satellite*'s true elevation angle 
and its apparent elevation angle (the arrival angle) as viewed from the ground 
station. Equation (26) is the equation for the elevation angle error; the rest of 
this report will be concerned with solving this equation, and with finding prop- 
erties of the function 8 (a. , R ) implicitly defined by it. 

g t S 

It should be stressed that all the operations indicated in Equation (26) can be 
carried out, at least in theory. From our previous analysis we know that, for 
any e -fimction in the class £ * and for any point (a^ , R^ ) e D, the following 
thills are true: 

• [1 - (1 + £ )^ (R^ /R)^ cos^ a^3 ^ 0 for all e [0, tt/2] and for all R ^ R^, 
with equality holding only for = 0 and R = R^ . 
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• The integral on the left-hand side of Equation (26) exists for all e [0, n/2] . 

• Equation (26) has a unique solution for , which lies in the interval 
[0, 7t/2); equivalently, it has a unique solution for , which lies in the 
interval [-a^, 77 / 2 -aj. 

We therefore do not need to worry about whether the square root is real, or 
whether the integral converges, and we know that the function ^ 

fined at every point of the region D in the plane. Practical difficulties 

in computation will arise, however, for values of close to zero, because of the 
fact that the integrand (on the left-hand side of Equation (26)) becomes infinite 
as R-* R^ , for = 0. (For realistic refractivity profiles will be very small, 
so that will differ only slightly from .) 

It will be convenient to write Equation (26) in an abbreviated notation. To 
this end we define the functions 




by 




Re 

(1 + e) —cos (a + /3) dR 
R2 


/Re 

v2 I 


' - (1 + ey j cos2 (a^ + /3) 


1/2 


for each e e £ ^and each p eD; the g-functions corresponding to e (R) = 0 will be 
denoted by go.p. (The symbol £ denotes the real number system.) With this 
notation we can write Equation (26) as 


g*,p(Sg) = go.pC®)- 


(27) 
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We note that: 


1* g«,p (/3) decreases (strictly) mono tonic ally to the value zero at /5 -tt/ 2 • 
; since < tt /2, it follows that ^ (0) > 0. 

2. If e, (R) > £ 2 ^^) equality not always holding, then g^ (/ 3 ) 

^ g, (/ 3 ) for all S e [-a^, tt/2 - o-^) . ^ 

3. g*,p (0) - go, p(0), with equality holding only for e (R) ^ 0. 


All three of these properties hold for any e -function in 8 ^ and for any point p 
in D. The third property can either be derived from the second, if we note that 
£ > 0 for ail R, or it can be deduced from the first property and the fact (which 
we already know) that Equation (27) has a unique solution. 


Some Properties of the Solution 


Positivity - We are now ready to derive some properties of the function 
8^(a^ , R^). We already know that this function is uniquely defined on the region 
D. We will now show that it is everywhere positive: for every € -function in the 
class 8 except the function £ (R) - 0, 


0 < 8 < — _ 
B 2 


(28) 


for allpeD. (Clearly, for = 77 / 2 , 8 ^ = 0.) This follows from properties 1 
and 3 of the g- functions; it is obvious from Figure 4. Since we now know that 
8^ > 0, we can limit the domain of definition of the g-functions to the interval 
[0,77/2 - a J . 

An Upper Bound -- The quantity ( 77/2 - ) in (28) is an upper bound for the 

-function; it is a very crude bound, however, and we will now derive a sharper 
one. If we ’’solve” Equation (26) for the 8^ in the numerator and delete the 8^ 
inside the square root, it follows from the positivity of 8^ thi’i 
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for any e e 8 ^ which is not identically zero and for all p eD. (Note that the ratio 
of the two integrals is less than unity, since e (R) > 0, and therefore Bj > 0.) Un- 
fortunately the bound is not too useful, because it requires the computation of 
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an integral and because it is not very sharp, although it is much sharper than the 
bound (tt/2 ). A computer program (which will be described in a subsequent 

report) was used to evaluate the bound for the fairly realistic refractivity 
profile described in the appendix; this profile will be referred to as the ’’Sperry 
profile," because it was taken from a report of the Sperry Gyroscope Company. 
The results Indicate that is quite sharp for large values of a, , but becomes 

progressively worse as decreases; for small values of it is very crude. 

For the Sperry profile, and for a signal frequency of 140 megahertz, B^ exceeds 
the true value of Sg by only a few percent for > 75°; for = 40° the er- 
ror is already about . 00%, and for = 10° B^ is about ten times larger than the 
true value of . (T>.e calculations were done for a satellite height of 1,000 
kilometers.) An attempt was made to make Bj more useful by replacing e (R) 
by a simple majorant function consisting of straight-line segments (of the form 
given by Equation (25)) and explicitly evaluating the integral. The attempt failed 
on two counts: the resulting formula was undesirably complex (although the 
majorant function consisted only of two line segments, an initial sloped line fol- 
lowed by a horizontal line), and the sharpness of the bound was severely de- 
graded, the new upper bound being almost 200% in error even for = 87°. For 
the Sperry profile, which is a fairly realistic one, the true value of varies 
from 3.9 milliradians at = 5° to 0.026 milliradians at 87°, for a satellite 
height of 1,000 kilometers and a signal frequency of 140 megahertz. As in all 
the other calculations made with the computer program, it was found that 8^ 
decreases monotonically with . 

Another Upper Bound - A different upper bound for can be easily de- 
duced from Equation (26); it is rather crude, but it has the advantage of great 
simplicity. To find it we just note that, formally, the left-hand integrand in 
Equation (26) differs from the right-hand integrand in only one way: the appear- 
ance of the quantity Q = (1 + €) cos (a^ + S^) in place of cos . Since the left- 
hand integrand increases with increasing Q, for any value of R, it is clear that 
Equation (26) cannot possibly hold if the minimum value of Q is greater than cos 
or the maximum value of Q is less than cos , In other words we must have 


and 


(l+mine^)-cos(a^+S^) < cos 


( 1 + max ) • cos (o-t + Sg) > cos 
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for any t -function which is not identically zero on the interval [R^ , R J , where 
min 6^ and max e^are the infimum and supremum, respectively, ofe (R) 
on the interval [R^ Since min = 0, the first inequality just implies 

that ; ^ > 0, which we already know. If we solve the second inequality for 5^ , 
taking into accoimt the fact that + g < tt /2 (see Equation (28)), we obtain the 
upper bound 


S < cos 

g 


-1 


cos a. 


1 + max e 


- a 


t ’ 


where the arc cosine is taken in the first quadrant. This bound holds for any 
point (a^ , Rj ) e D, and for any e e which is not identically zero on the interval 
[R^ , R^] . We can get rid of the dependence on R^ , at the cost of weakening the 
bound slightly, by replacing max e ^ by the possibly larger number , the su- 
premum of € (R) on the interval [R^ ,oo). We then have 


0 < S < cos~^ 

K 


cos \ 


(30) 


= ^2 (a,, R,), 


for any e -function in £ ^ which is not identically zero on [R ^ ) , and for any 

point p£ D. 

The bound B 2 was tested numerically for the Sperry profile. The results 
indicate that the ratio of to the true value of is remarkably constant over 
a large range of elevation aisles; for the profile used, and for a signal frequency 
of 140 megahertz, this ratio remained between 1.79 and 1.83 for between 20® 
and 87°, for a satellite height of 1,000 kilometers. For < 20° the ratio in- 
creased; at 5° it was 2.75. It would appear from these results that B^ is useful 
as an easily computed estimate of 5^ . Further, the relative constancy of the 
ratio Bj /S indicates that might be used to form approximate formulas for 
S by multiplying B^ by a polynomial of low degree in (and possibly R^ also); 
such a formula might be very useful over a given range of elevation angles, for 
a specific refractivity profile. For the profile used in the tests, it would be suffi- 
cient to divide B 2 by 1.81 to obtain a formula accurate to about 1% over the entire 
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Table 1 

Computed Values of Elevation Angle Error, in Milliradians 


Ele- 

vation 

Angle 

(deg) 

Height 

(km) 

100 

200 

300 

400 

500 

1,000 

5,000 

10,000 


1,000,000 

1 






8.00 





2 






6.33 





3 






5.23 





4 






4.46 





5 





3.89 

3.89 


3.84 


3.83 

10 






2.42 





15 





1.88 

1.74 


1.55 


1,51 

20 






1.33 





25 





1.21 

1.06 


0.863 


0.830 

30 






0.856 





35 





0.845 

0.706 


0.548 


0.524 

40 






0.588 





45 

0.290 

0.301 

0.435 

0.599 

0.604 

0.492 

0.385 

0.370 

0.356 

0.354 

50 






0.412 





55 





0.428 

0.343 


0.253 


0.242 

60 






0.283 





65 





0.287 

0.228 


0.166 


0.159 

70 






0.178 





75 





0.165 

0.131 


0.0947 


0.0906 

80 






0.0860 





85 





0.0540 

0.0427 


0.0308 


0.0295 

86 






0.0341 





87 






0.0256 





88 






0.0170 





89 






0.0085 






range = 20° to 87°, for the satellite height of 1,000 kilometers and the given 
signal frequency. It would not be difficult to account for the variation of 5g with 
Rji , since Sg varies much more strongly with than with . To illustrate 
this, table 1 shows the values of (obtained with the computer program) for a 
range of values of elevation angle (a^) and height (R^ ~ Sperry re- 

fr activity profile and a signal frequency of 140 megahertz. Some of the data in 
table 1 are plotted in Figure 5. (Many of the positions in the table are blank be- 
cause those values of were not computed in the test runs.) 
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5-El.votio„ Anal. E^o, as a Function of Elovotion Anglo, fo. Vorious Altitud.s 


A Lower Bound - So far we have shown that S > 0 and we have derived 

g 

upper bounds for It, It Is of some interest to also derive a lower bound, other 
than zero. To accomplish this we start again with Equation (26), and rewrite it 
in the form 





Re ^ 
R2 " 

os dR 

/R 



r 2 ~i 

COS'* 


1 c ^ /*T 


t 1 

\R 

1 U a 

T 

_(1 + e)2 cos^ (a^ + 8^) 


1/2 



'2 


Since (a^ + 8 ) < tt/ 2 and < n 1'^ for all points PeD, we see that neither of the 
cosines is zero; the above equation, whose left-hand integrand was derived from 
that of Equation (26) mainly by multiplying and dividing by the two cosines, is 
therefore equivalent to Equation (26). We now apply the mean-value theorem to 
the function to obtain 



fKV 1 

f * \ 2 


cos^ 

\ 1 cos Or. 

\ R / 

+ 

_(1 +6)2 cos^ (a^ + Sg) 


1 


1/2 




where 0 < 5 < 1 and 6 depends on R. This is valid only for > 0; if = 0, 
the quantity [1 - (R^ /R)^ cos^ ] becomes zero for R = R^ . If we substitute 
this into our rewTitten equation and simplify, we are left with the equation 


35 



cos a. 


f 




cos r 




R \ 


j COS^ 


r. 

[ (1 ■! cos^ (a^ + ^ ) 

■’T 


- 1 


dR 


= 0 . 


cos“ a. 


(1 + t )2 cos^ (a^ + 6^) J 


- 1 


3/2 


We will denote the bracketed quantity in the numerator by Q, for convenience. 
From the relation (30) it follows that 


cos^ 

(1 + cos^ (a^ + 5^) 


^ 

(1 + 6)2 


Since the right-hand side of this inequality is unity for any value of R for which 
c taltes its maximum value, we see that the left-hand side is sometimes less 
than unity, whence Q < 0. On the other hand, for any value of R for which 6-0, 
for example for R = R^ , the left-hand side is clearly greater than unity (since 
Sg > 0), whence Q > 0. In other words, Q is sometimes positive and sometimes 
negative. Now suppose we set 6?(R) = 0; clearly, the Integrand will then become 
larger in magnitude for those values of R for which Q is positive, and it will be- 
come smaller in magnitude for those values of R for which Q is negative. It 
follows that 


f 


K 

R2 


cos* a. 


• cos 


- 1 


dR 


(1+ 6)2 cos2 (a^ + S^) 


> 0 


1 - 

L \R 


3/2 


cos* 


If we solve this inequality for we get our lower bound: 


cos \ 

/TTa j 


- a. 


r.L(a^.R,)>0. 


( 31 ) 
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where 


f 


A r: 


dR 

R2 

/r y 1 

1 - — 1 cos^ <1 

\.R / ‘ 

3 2 


- 1 


I 


dR 


** R^ (1 + e)2 I 1 J 


3 '2 


It is clear that A > 0, and therefore L > 0. This bound holds for any e -function 
in 8 which is not identically zero on the interval [R^ » 1^ s ] for any point 
pe D for which > 0. It was tested numerically with the computer program; 
the results indicate that it is quite accurate. For the Sperry profile, at a signal 
frequency of 140 megahertz and a satellite height of 1,000 kilometers, the lower 
bound L differs from the true value of by only about 0.7% at =5°. For 
larger values of the error steadily decreases. At = 15® it is 0.2%, at 25® 
only 0.1%, and for > 45® it is less than 0.067o. It thus would appear that the 
lower bound L, as given by Equation (31), is a rather accurate approximation for 
the elevation angle error 8^ . 


Partial Derivatives - Co far we have proved that the solution function 
5^ (a^ , ) is uniquely defined on the region 


D - R.) I 0 < < n/2 and R, > R^ } - [0, tt/2) x (R^ , oo) 


of the - R^ plane, for any e -function in the class 8 , and we have derived 
some upper and lower bounds for it. We will now try to get more detailed in- 
formation about the behavior of the solution function. Some information can be 
obtained with the help of the implicit function theorem of analysis. Suppose we 
write Equation (26) in the form 


F(Sg; s-R,) = o 


bv iitting all terms on the left-hand side of the equation. The function F, which 
maps a region E of the 3-dimensional space into the reals, is 

given by 
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(1 + e) —cos (a + & ) • dR 

R2 * 




1 "' 


r 


— cos a • dR 
R2 


-4 


cos' G. 


1/2 


for ( , R^ ) g E. The region E is given by; 

E - {(S^, , R^) e 1 e [0, n/l) and R^ e (R^ , '^) and e [0, tt/2 - a^)} 


For easy visualization this region is shown in Figure 6; it is a semi-infinite 
prism, with a cross-section in the form of an isosceles right triai^le. The 



Figure 6-The Region E 
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region E does not include the entire surface of the prism; it includes only the 
interiors of the vertical and horizontal faces plus that part of the -axis for 
which R > R^ • The base of the prism (i.e., the horizontal face) is just the 
region D. We will denote the interiors of E and D by E and D, respectively. It 
is clear from our previous results that the graph of the solution function 
S » R s ) is a surface which lies in E for any e -function in which is not 
identically zero. (If e = 0, the surface just becomes the region D itself, since 
^ 0 for all points of D in this case.) We could have defined the function F at 
all points of the slab [0, tt/2) x (R^ ,(») x (_oo, oo); but then F would be periodic in 
h . Since we already know that the solution surface lies in the region E, there 
is no point in complicating the problem by defining F over a larger domain. 


We are now in a position to apply the implicit function theorem; we will 
state this theorem as it applies to our particular case. Suppose the following 
hypotheses are satisfied: 


1. The partial derivatives BF/Ba^., BF/BR^ exist and are continu- 

ous on the domain E. 

2. There exists a point (8^, , Rq )e E for which F (8^, R^) = 0. 


3. 


BF 

BS 


^ 0 . 




The theorem then states that there exists an open neighborhood W C D of the 

O 

point (a^, Rq) e D and a unique function h; W - R, such that: 

1. h(a^.R,) = 8 q. 

2. F(’ a^,R^)=0 for all points (a^ , R^) e W. 


3. 


—11 , exist and are continuous on W. 

Ba/ BR 

t s 


BF 


Bh 

Ba^ 

(h(a^.R,),a^.R,) 

Ba 

t 

(a . R, ) 111 

BS 

K 

(h{a,.R,). a^.R.) 
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and 


3h 


. IL 

^^8 (h(g^.R, ),a^.R, ) 

"aF 

(h(a^.Rj, a^.R^) 


for all points (a^ , ) e W. 

In other words the theorem states that, if we know a point (Sq » which 
the equation F (8^ , ) = 0 is satisfied, then this equation can be uniquely 

solved for in terms of and R^ for all points (a^, R^ ) sufficiently close to 
(ttQ, Rq ). It is important to note that the theorem guarantees the existence of a 
continuously differentiable and unique solution function only locally; even if hypo- 
theses 2 and 3 are satisfied for every point in f), it does not follow that there 
exists a global solution unique on D, although there exist locally unique solutions 
about each point of D. However, we know from our previous work that, for each 
point Rg )e D, there exists a unique such that F (S^ ,a^ » With 

this additional information we can conclude that the solution function h, or S 

O ^ 

(a^. , R^ ), is defined uniquely everywhere on D and has continuous partial deriva- 
tives on D. It remains only to show that hypotheses 1 and 3 are satisfied. It is 
easily seen that the three partial derivatives exist and are continuous on E; they 
are given by the following formulas: 


(1 + e)-!-dR 
R2 


1-(1+€)M — ) cos^ (a^ + sp 


dF . , 

- — = - sin (a 
da^ '■ 




R 

(1 + e)-ldR 

R^ 

/r \2 -|3/2 

1 - (1 + e)2 j cos2 (a^ + 5^) | 
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BF 

BR 


[1 + e (R )] —cos (a + h ) 
R? " 


R 


R2 

S 


COS 


1 - [1 + £(R^)]2 



— 

1/ 2 

1 

(N 

od 

1 

1 2 

cos^ (a^ + S ) 

' t ^ g/ i 


1 - j cos 


_J 


_ \ s / _ 



It is clear from the formula for B F/BS^ that this derivative is non-zero (in fact, 
it is negative) at all points of We are therefore assured that the solution 
function Rj), which we already know exists and is unique, has continuous 

partial derivatives on 6; by conclusion 4 of the theorem, these partials are given 
by the following formulas: 


dR 


B8 

i 

Ba. 


sin a 

“ 1 1 • 

' R2 

' fK y 1 

1 - ( — ) cos^ a 
\R / J 

3/2 

' sin (a^ +8g) .R, 


(1 + e) dR 



(32) 


R2 


R \ ^ 

1 - (1 + ( Y ) (""t + ^g> 


3/2 


R, 

[1+ € (R )] — cos (a +2 ) 

R2 ‘ 


COS a. 


BS 

8 

BR. 


1 - [1 + 6(RJ] 


/R y 
^ (r” ) 


1/2 




2 1 

1/2 

1 - 

e 

cos^ a. 







sin (a 




(1 + e)_ldR 
R2 


' 1^1 - (1 + e)^ j cos^ (a^ +8g) J 


3/2 


(33) 

In these formulas, the value of must be taken as (a^, R^), the actual solu- 
tion value of 8 at the point ( R^), in order to evaluate the partials at the 

point (a , R ).* Although we know that the solution for 8 exists and is unique, 
we do not know what it is. In this sense the above formulas are incomplete; to 
render them complete we would have to kiiow the solution function. In spite of 
this, we can still get some information from these formulas, if we use what we 
know about the solution function. 
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It should be pointed out that much more information is obtainable from the 
implicit function theorem than we have stated above. For example, it is easy to 
see that the function F is differentiable any number of times; i.e., the partial 
derivatives of all orders exist and are continuous. The theorem (in its more 
complete form) then allows us to conclude that the solution function (®t » ^s) 
has continuous partial derivatives of ail orders. We will only use the first-order 
paiiictls, however, because enough information is obtainable from them to get a 
general idea of the behavior of the solution function and because the higher-order 
partials are given by excessively complicated formulas. 

We will examine the partial /BR^ first, since it turns out to be the easier 
of the two to handle. If we look closely at the two terms in the numerator, we 
notice that the first term differs from the second in only one way: the quantity 
Q = [1 + e (Rg ) ] • cos (ttt + 5g) appears in place of cos ; furthermore, the 
first term increases monotonically with Q. It follows that the algebraic sign of 
is just the sign of Q - cos o .^ ; i.e., 

sign (BSg 'BR^) = si gn { [1 + e (Rg)] ‘ cos [a^ + (a^ , R^)] - cos } 

^sign{k (R)>. 

t ^ 

From the upper bound for given by the inequality (30) it follows that 

(1 + ~ ^ 

from the definition of as the supremum of e(R) on the interval [ ,qo) and 
from the continuity of e and 8g as functions of R^ it follows that, for any given 
value of , there exists at least one point in (R^ ,oo), and an open interval con- 
taining that point, such that the function k(R^ ) is positive in that interval. In 
other words, 

BS 

— ^> 0 for R z (34) 

BR ® 

S 

where 0.* is some open set. contained in the interval (R^ ,®). (Note that the set 
fi'*' depends on the angle .) We will now show that (R^ »°“). To see this, 
we just recall thate (R^ ) = 0; since > 0, it follows that k^^ (R„ ) < 0 for all 
. It follows from the continuity of the k-function that k^^ (R^ ) is negative in 
some open neighborhood of the point R^ , and possibly on other open intervals 
as well. In other words. 


BS 

0 for R c Q', (35) 
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where is some open set containing a neighborhood of the point and con- 
tained in (Rg ,00 ); note that Cl~ depends on . From all this we see that, for 
each e(0,7r/2), the partial derivative is sometimes positive and 

sometimes negative, and therefore also sometimes zero. This behavior is 
evident In the data shown in table 1; the table also indicates another property, 
which is clear from Equation (33): 




BR. 


0 as R 


(36) 


for any e ( 0 , 77 / 2 ). 

We now turn to the partial /Ba^. From values of obtained with the 
computer program it is seen that, for the Sperry profile, decreases monotoni- 
cally with for any given value of R^ (e.g., see Figure 5). Since there is noth- 
ing very special about the Sperry profile, it seems reasonable to suppose that 
this monotonicity might hold more generally. We will investigate this question. 

It is evident from Equation (32) that a sufficient not necessary) condition 
for the derivative to be negative is that the ratio of the two integrals be less than 
unity: 



(1 ^ €) dR 

/R \2 “ 13/2 

l-(l + e)2 cos2(a, 



dR 



[ fR y 1 

3/2 

R2 

1 - [ — ] cos^ a 



L \ R / J 




(37) 




(Tho ntegral was evaluated by using the substitution (R^ /R) cos = cos &.) We 
can get rid of the unknown S and at the same time construct conditions which 

g 

imply (37) by replacii^ 5 in (37) by any of its upper bounds. Our crudest and 
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simplest upper bound for is the one given by (28): < v /2 - a.^ . The cor- 

responding condition is 



It is evident that $ is a monotonically decreasing function of a-^ . Further, O(a^) 
becomes infinite as - 0, and 






% 

(l + e)dR 

R 

e 


We conclude that, for any e -function in S ^ , there exists a value of for each 
value of , say a" (R^ ), such that BS /Ba.^ < 0 in the interval (a" ( ), 'n/2) 

(and probably in a larger interval), w!ierea"(R^) is the value of for which (38) 
becomes an equality. 


The value of a" given by Equation (38) would probably be close to 90®; we 
can get a smaller value of oT by using a sharper upper bound for 8^. If we use 
the bound B ^ given by Equation (30) , we get the condition 



1 


(1 +e) dR 


(1 t 

(1 + € )2 



COS* a. 


3/2 




(39) 


We reach the same conclusion as before concerning the existence of a* , because 
the left side of (39), evaluated at = tt/ 2, is clearly greater than <l>( 77 / 2 ). The 
value of a" computed from (39) will in general be smaller than that computed 
from (38). 

More detailed results can be obtained, but they are probably not worth the 
considerable effort required to find them. For example, after a very lengthy 
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calculation (which will not be given here), the writer found that a“ (R, ) - 0 as 
R^ - 00 ; however, all attempts to prove that a" (R^ ) 0 (i,e., that the derivative 

d? /da, < 0 for all a, and ail R ) failed. 

t t S 

We will now go on to a different topic - the transformation of Equation (26) to 
a form which is more convenient for numerical calculation. 


A Transformed Version of the Elevation Angle Error Equation 


Earlier we wrote Equation (26) in the abbreviated form (27), by using the 
functions ^ (/3). It would be useful to know something about the derivatives 
of these functions. We find by direct calculation that 


gf.p (^) ^ - K ■ 

and 

(/5) = - Re ■ 


f 


(1 + e) dR 




1 - (1 + e)2 ( J cos^ (a^ X (3) 


r 


(1 + £ ) dR 


3/2 


R 

R2 !l-(l + e)^l— I cos^ (a^ + /3) 
L V R 


3/2 


+ 3R^ cos (a^ + f3) • sin ^ (3-t + /3) 


r 


(1*0* 1 ^) dR 




1 - (1 ^ I y) 


[5/2 


The first derivative is clearly negative, but the behavior of the second derivative 
is not very clear. Higher derivatives will be even more complicated and diffi- 
cult to analyze. It seems that the functions (/3) show undesirably complex 
behavior. Fortunately, it is possible to avoid tfiis difficulty by making an appropri- 
ate transformation of the dependent variable in Equation (26). To this end we de- 
fine the functions 


cp^ : [0. n/2 - aj - [1, <») 

t 
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by 


, cos a. 

cp ; h " ^ . 

cos'* (a^ + p) 


for any value of in the interval [0, tt/2). In terms of the variable i , we can 
define the functions 


by 






(1 + e)_ldR 


R 

^ - (1 + €)2 \-^ I S 


1/2 


It is convenient also to change the variable of integration to p , where 


^ R 


(40) 


If we define the function e {p) by 


e (/5) = e (R) = € (R ■>) , 


(41) 
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we can write the f-functions in the form 


^,p( 0 = J 

Jo 


dp 


I ^ - p^cos^ a 

( 1.^)2 


12 


(42) 


where =R^/Rjj . Clearly e (0. 1). 

We can now write Equation (26) in the form 


f 


dp 


i cos^ 


.(1.6)2 


1 / 2 




dp 


p2 cos2 


cos"^ (p^ cos a^) - 


cos a. 


(43) 


or 




where f is the f- function corresponding to e (R) = 0 and 


cos2 

^ 

cos2 (a^ + 


We also write the inverse relation for reference; 



(44) 


(45) 


( 46 ) 
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To avoid confusion, it should be noted that Equation (45) defines c, not as a 
function of , but as a function of and (since depends on and ). 
However, from the fact that is a monotonically increasing function, it is 
clear that any upper (or lower) bound for generates an upper (or lower) bound 
for ^ , and vice versa. For example, the bounds for ^ , corresponding to the 
bounds for given by (30), are 


1 <<i ■ 


( 47 ) 


These bounds can also be read off from Equation (43) by inspection, if we recall 
that? (p) > 0 for all p . (Note that this upper boimd does not depend on .) 

From (47) we conclude that the quantity ^/[l + e(p)]^ is less than unity 
for some values of p and greater than unity for other values oi p. To see this 
we note first that, for any value of p for which ?(p) = 0 (e.g., for P = 1), 


[1 


I 


£ 

1 


> 1 . 


Secondly, we note that for any value of p for which ? (p) = * 


[1 



(1 + 


< 1 . 


Of course it is also evident from inspection of Equation (43) that this must be 
so; if the above quantity were either always less than or always greater than 
unity Equation (43) could not be satisfied. 


Equation (43) is our new form for the elevation angle error equation. We 
will now show that it has more convenient properties than does Eqxiation (26). If 
we look at the derivatives of the f-functions, we easily find that the n-th deriva- 
tive is given by the simple formula 


= (-ir 


(48) 


1 • 3 • • • • 

•(2n-l) r 

dp 


2" 


^ 2 2 ' 
^ p* COK* a. 

n + l/2 
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for n = 1, 2, 3, This very clear formula contrasts sharply with the ob- 

scure behavior of the derivatives of the g-functions. In particular, we note that 
the first derivative is always negative and the second derivative is always posi- 
tive. We will use this fact later, when we describe how to solve equation (43) 
numerically by the Ne^vton-Raphson method. 


If we denote the integral on the right-hand aide of Equation (48) by 1(0. 
where n = 0, 1, 2, 3, • • * , it is easy to see that for any ^ -function in ^ * : 




— - cos* a 

all C [1, 3o)and all p < 
C 1 and p = 1; 


^ 0 for all e [0, tt 2 ) , 

1, with equality holding only for = 0 , 


• p(C ) (" exists (i.e., is finite) for all [ 1, ou)^for any point p^^ D; 

• I„(0» and therefore also f (O, exists (i.e., is finite) for all C [1. 
and for n = 1, 2, 3, • • • , for any point pe D; 

• I„(l), and therefore also f (1), where n = i, 2, 3, * * * , does not exist 
(i.e., is infinite) for = 0. 


The first three properties follow directly from our previous results (see the dis- 
cussion following Equation (26)). To see that the fourth property holds, we 
proceed as follows: 


In(l) 


=r 


dp 




2n 


1 


2 2 
-p cos^ 


n + 1 ^2 


J -1 

I P dp 

[1 - p* cos* ] * ■*■* '* 


(1 + o 


*" cos* 


sin a. 


/T 


1 

p* COS* 
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Clearly, the last expression becomes infinite as - 0. The fact that the 
derivatives become Infinite (at C = 1) as - 0 causes difficulty in the numeri- 
cal solution of Equation (43), for small values of . 

Tt is clear that Equation (43) is equivalent to Equation (26), for all € -functions 
in t and for all points peD; thus the solution of one of these equations can be 
obtained directly from the solution of the other by using Equation (45) or Equa- 
tion (46). From this point on we will use either Equation (43) or Equation (26), 
whichever happens to be convenient at the time. 


A Method for Solving the Transformed Equation Numerically 


In terms of the f-functlons (defined by Equation (42)), the transformed ele- 
vation angle error equation is given by Equation (44); we will write that equation 
in the simplified form 




(49) 


for convenience. A plot of this equation will have the general appearance shown 
in Fig^ure 7. 
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Of the various methods available for solving such an equation, we will 
describe only one: the Newton-Raphson method. The reason for this choice is 
the fact that f (O ^ 0 and f* (O > 0, for all ^ e [ 1, oo); this guarantees that the 
Newton-Raphson process will converge monotonically upwards to the solution 
value provided that the starting point is smaller than ^ . We will now describe 
the process in detail. 

The Newton-Raphson method uses a succession of approximating tangents to 
the curve y = f(C). The first step is to pass a tangent line through the point 
where ^ ; the equation of this line is 


y = (^o) • (^ -^o) + f (^o)- 


The value of C for which this line crosses the horizontal line y = fg gives us the 
first iterate, ; .ve calculate by setting the above expression e<^al to f ^ 
and solving for C • The result is 






f (^o) - f. 

-f' (f„) 


To get the second iterate, * we pass a tai^ent line through the point , f (^ j)) 
and calculate its intersection with the line y = fg ; we find that 


^2 ~ 


f tfi) - h 


Clearly we have the general recursive relation 






(50) 


for n = 1, 2, 3, * • • . 

To start the computation we need a value for^g which is less than the true 
solution The simplest choice is to take ^g =1; this will always work in 
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principle, except whena^ =0. In practice, computational difficulties will arise 
even when is just close to zero, because of the boundless growth (as - 0) 
of f '(1) and of the integrand of f(l) at p - 1. Somewhat better results can be 
obtained by taking to be a lower bound for other than unity. For example 
the computer program (which will be described in a companion report) takes 
to be the lower bound for f corresponding to the lower bound for S given by 
(31). 


The iterative process is stopped when some convergence criterion is satis- 
fied. Instead of putting a condition on the iterates , the computer program puts 
a condition directly on the corresponding values of » which we will denote hy , 

because is the quantity we are really interested in. The criterion used is 
that the relative error ( 8^ + j )/S^ ^ , be less than a preset value; when this 
condition is first satisfied the iterative process stops. 

In general, the Newton-Raphson method converges rapidly after the iterates 
get close to the true solution. But during the first few iterations some other 
process, such as the secant method, may be more efficient, and it may be ad- 
vantageous to use such a method to start the computation. In our case, however, 
the iterates converge monotonically upwards toward tie solution, and it seem*! 
simpler to start right off with the Newton-Raphson method (which is what is done 
in the computer program). 

To complete this discussion of the numerical solution of Equation (49) we will 
prove by induction that the iterates do in fact converge monotonically upward 
to the solution , and we will show that a sufficient condition for this to happen 
is that f‘ii)< 0 and f " ( O ^ 0 for all C . 

Suppose then that ^ ^ ^ given by Equation (50). If we 

could prove that 


^ n - 1 




(51) 


then we could conclude that the sequence (^ ) is monotonically increasing and is 
bounded above by ^ . (The induction is started by noting that < ^ , so that the 
premise is in fact true for n = 1.) To prove the first inequality in (51) we use 
the fact that f(C) is monotonically decreasing. From our premise that ^ ^ 
it then follows that f(^ri- 1 ) ^ f(^) ” ^ • With this result we can conclude directly 
from Equation (50) that i , if we also use the fact that f'(^„_ j) < 0. To 

prove the second inequality in (51)- we use Taylor's Theorem to write 
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where - 1 < ^ . If we use Equation (t'iO) to replace the first two terms on 

the right by f^ we obtain the relation 




It follows that f(^^ ) > f^ = f(^), because f *(^*) > 0. From the monotonicity of 
f(0 we conclude that which proves (51). 

There is one last thing to prove - we must show that ^ . We have just 
shown that the sequence (^^) is monotonically increasing and is bounded above 
by i ; the sequence must therefore converge to some number 5 . To show 

that = i we observe that, since is convergent, 0. From 

Equation (50) we then conclude that f(f^) - f^ . Since it is clear that the inverse 
function f ~ ^ (which exists by virtue of the monotonicity of f) is continuous, we 
conclude that i . 


Some General Formulas Giving Approximate Solutions to 
the Elevation Angle Error Equation 


General Form — We approach the problem of finding approximate solutions 
in the following way. First we obtain, by any means available, an approximation 
to the solution of Equation (43) in the form 


^ = 1 +A, 


(52) 


where 0 < A « l. The fact that A is small follows from the bound given by (47) 
for i , namely ^ < (1 + e^)^ ; since will be small (of the order of 5 x 10“^) for 
a realistic refr activity profile, we see that 


A = ^-1<2€ 

^ m m 
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Nexv we write the correspondii^ approximation for by using Equation ( 46 ); 


b - COS 

e. 


cos a. 


/TTa 


- a 


i ■ 


( 53 ) 


If is not close to zero we can simplify this formula by using a differential 
approximation. We regard the expression on the right as a function of A, say 
A(A), and consider as a parameter, ignoring the fact that A actually depends 
on . By using Taylor’s theorem we can write 

A(A) = A( 0 ) + V (t?) • A 


— cos a. 
2 ’ 


I 


A, 


cos‘ 

1 . (1 + 7j)3/2 


1+77 


where 0 < 77 < A. If is not too small we can neglect 77 without affecting the 
last expression too much, because 77 is small; in that case we may write 


S 

g 



cos a 

i • A. 

s in 


( 54 ) 


We can get some feeling for the range of validity of this approximation by writing 


cos'^ a 


1 = 1 - (1 - 77) CO a = s in^ a + 77 cos ^ a 

1+77 


so that 


1 cos a. 


2 s in a. 


/T 7 


77 cot^ 
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The neglect of rj will be justified if 17 cot^ « 1 ; in this case we can write 


A 



cos 

sin 


1 

2 ^ CO 


t2 



• A. 


From this expression we see that Equation (54) will be accurate to about 10% or 
better if we have 1/2 rj cot^ <0.1. From the comments following (52) we see 
that this last condition will be satisfied if e cot^ < 0.1. Since e = 5 x lO"^ 
for a realistic refractivity profile, the condition becomes tan > 1 / 5 * 10 "^, or 
roughly > 12®. In other words, the value of given by Equation (54) will 
differ from that given by Equation (53) by no more than 10% or so, for > 12®. 

Before we go on to obtain explicit expressions for A , some comments should 
be made about the problem of determining the accuracy of the corresponding 
formulas for 8 ^ . When the various formulas were first derived, error terms 
were included and attempts were made to use them to evaluate the accuracy of 
the formulas. This approach turned out to be intractable. The trouble is that 
the error terms are quite complex, involving complicated definite integrals. 
Attempts to simplify them by using estimates or bounds invariably result in ex- 
pressions which give far too crude an estimate of the error. The original error 
terms could of course be evaluated numerically, but it is simpler and more 
direct just to compare the approximate formulas with accurate numerical solu- 
tions of the elevation angle error equation, and this is the course we will follow. 

First Newton-Raphson Iterate - Perhaps the most obvious candidate for an 
approximate general formula is the first iterate in the Newton-Raphson numeri- 
cal solution scheme which was described earlier. If we take = 1, it is seen 
from Equation (50) that the first iterate is given by 




= 1 


+ —————— , 
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in other words, 


A = 


f 

-f' (1) 


f ^ f' 

Jp^ [1 - (1 + cos^ aj [1 _ 


dp 
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P ^ cos^ 


i f' 

^ •'■». (1 t e)^r !- 


dp 


L (1 - O- 


2 2 
- cos^ 


3/2 


It is not practicable to go further than the first iterate; higher iterates would in- 
volve integrals within integrals. 

A Differential Approximation - Aside from using the first Newton-Raphson 
iterate (or the first iterates in other numerical schemes), the most obvious way 
of obtaining approximate formulas is to use differential approximations. To 
proceed in a systematic way, we use the notation of Equation (49): 


f (a = 


i 


dp 

^ - P^ cos^ a 


1/2 


(56) 


y\2 


(1^ e) 


We also consider the integrand as a function of ^ : 


r(^) = 


P 2 2 

- p cos^ 


1 /2 


(57) 


.(1.6)2 
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The elevation angle error equation can then be written as 


f (O = 



r (i) dp - . 


(58) 


The various possibilities for differential approximations are now apparent. 
The most obvious are the first-order Taylor expansions of f(^) and r(^) about 
the point ^ = 1 : 


f (0 = f ( 1 ) + f'(l) • (^- 1 ) 
and 

r (0 = r(l) + r'(l) • (^- 1 ). 

It is easy to see that if either of these expressions is put into Equation (58) and 
the resulting equation is solved for f , we just get back the first Newton-Raphson 
iterate. (It is also geometrically obvious that this is tri'.e.) 

Another possibility is to write y /{I +^)^ , and to regard r as a function 
of y , ignoring the fact that y actually depends on p . Since we know (from the 
bounds for given by (47)) that y is close to unity, we can expand r to the first 
order in y about the point 7 1 ; it follows that 


r(0 


1 

[1-/0^ cos^ 


_£_-i 

1 . 1(1 

^ [1-/3^ cos^ 


If we put this expression for r(^) into Equation (58) and solve for A(= 1), we 

obtain 


/ 

jp 


dp 


[ 1 -/ 3 * cos* a ]^^* 

S ^ 


A = 


f- 

Jp ( 14 - 


- I 


dp 


(1 4 -e^* [1 - p* cos* a^] 


3/2 


(59) 
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It is interesting to note that this is precisely the lower bound given by (31). The 
smallness of A can be better displayed by writing the expression as a simple 
fraction: 


r‘ 2 

(1 [1 - cos^ 

7ZZIZIZZ 

Jp^ (I + 7)^ [1 - cos^ ^ 2 


Another Differential Approximation - We can obtain yet another approximate 
formula by regarding r as a function of the two variables and 7 , ignoring the 
fact that 7 depends on p. Since f is close to unity and 7 is small^ we can expand 
r to the first order in ^ and 7 about the point ^ = 1, ^ =0; it follows that 


r (f ) = 


f(s-l) 


[l-yo^cos^a^] 


3/2 


[1 - p- cos^ ) 


3/2 


This approximation should be fairly accurate (for all /d ) if is not close to 
zero (if = 0 and p - 1, the square root vanishes). If we put this expression 
for r(^) into Equation (58) and solve for A( = ^ - 1), we obtain 


A = 


r — 

*//5. [1 - / 


e d p 


cos^ ] 


3/2 


I 


dp 


^ [1 - p^ cos^ ] 


3 / 2 


(60) 


2 sin a. 


1 - 


p^ s xn 


L 


€ d o 


/i - p? 


^ [1 - p^ cos^ ] 


3/2 


p‘ cos- 
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where the integral was evaluated by using the substitution p cos = cos ^ . 

Note that this formula is not valid if is close to zero. In fact, for = 0 
formula (60) gives A = 0. We know that this is wrong, because it implies that 
^ = 1 for = 0, which is impossible by virtue of the inequality (47). 

A Simple Formula for the Elevation Angle Error - Any one of the three ex- 
pressions for A given by Equations (55'> , (59) or (60) could be used together with 
Equation (53) or Equation (54) to obtain an approximate formula for . It is 
evident, however, that the simplest of the three expressions for A is given by 
Equation (60). Further, the other two expressions for A can be reduced to it by 
using approximations based on the smallness of ?(p). For these reasons it 
would seem that this formula might be the best one to use. If we combine It 'vith 
Equation (54) the factors sin cancel each other and we obtain the follovving 
formula for : 


S 


g 


cos 
P, sin 

/l - Pj cos^ 



e d p 

[1 - cos^ 


(61) 


In spite of the fact that neither Equation (60) nor Equation (54) is valid for 
close to zero, numerical comparisons with the correct solution show that the 
formula given by (61) is remarkably accurate down to smsdl angles. (This is 
probably due in part to the fact that the sin factors cancelled out.) Table 2 
shows the percentage error of formula (61) for various elevation angles, for the 
Sperry profile and for an altitude of 1,000 kilometers and a signal frequency of 
140 megahertz. The % error is defined by 


% error 


value of computed f rom ( 61 ) - correct value of 
correct value of S 

g 


X 100. 


The data In table 2 show that formula (61) gives slightly high values for ^ , at 
3-jast for this particular case. The formula is seen to be accurate to witiun 1% 
do vn to elevation angles of less than 10°; even for = 1°, it is only 8% too 
high. This accuracy is considerably better than one might have expected from 
the way in which the formula was derived. 
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Table 2 

Accuracy of the Formula for 8 ^ 


Elevation Angle 
(degrees) 

% Error 

89 

0.05 

87 

0.05 

85 

0.05 

80 

0.05 

75 

0.05 

70 

0.05 

65 

0.05 

60 

0.05 

55 

0.05 

50 

0.05 

45 

0.05 

40 

0.05 

35 

0.05 

30 

0.05 

25 

0.06 

20 

0.09 

15 

0.2 

10 

0.4 

5 

2. 

4 

2. 

3 

3. 

2 

5. 

1 

8. 


Numerical tests were also run on formulas for S obtained using the 
other two formulas for A . The results show slightly greater accuracy, but not 
enough (in the writer^s opinion) to compensate for the greater complexity of 
those formulas. 

Before closing this section, we should show that tiie integral in formula (61) 
always exists. We use the fsmt that e (R) e 8^ , so that c (R) < (1 -a ) (R - R^)/ 
R^ for all R > R^, by the condition (25). In terms of tiie variidsle p = R^/R this 
condition becomes €(/?) < (1 -a) {l/p - 1), for all /? < 1. We use this condition 
to obtain a bound for the integral: 
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€ dp 


cos^ 


< (1 -cr) 


1 ( — - 1 dp 


Jp [1- 


213/2 


p"] 


for all e [0, tt/2 ] 


= (1 -a) 


r 


dp 


p(l (1 + p) 


k3/2 


< f' 

P, (1+P,)^''^ i. 


dp 


(1 - 

2 (1 - 0 -) (1 - ^ 


P, (1 +Ps) 


3/2 


< <». 


That is, the int^ral is finite. 


Some Comments on Methods of Finding Approximate 
Solutions for Particular Cases 


In the preceding section we were concerned with finding general formulas 
which are a[^roximate solutions of the elevation angle error equati(»i. In this 
section we will comment on the problem of finding simple approximate soluti(xis 
for specific refractivity profiles, with the aim of findii^ ways of minimizii^ the 
amount of computation needed to find S . Two possible approaches to this prob- 
lem will be described briefly. 

Numerical Approach - The basic idea is to solve the elevation angle error 
equation for various values of and (for some given signal frequency or 
frequency band) and then to fit curves to the computed points. Since varies 
less with R, than with , it would probably be feasible to obtain a number of 
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simple formulas for S in terms of , for various values of . It might also 
be feasible to fit a surface (i.e., a function of two variables) to the computed 
points. In either case, the aim would be to obtain simple, easily evaluated 
formulas for 5^. 

Analytic Approach — The idea here is to develop simple formulas for 
by using approximations for the integrals in general formulas such as formula 
(61). For example, one might use the trapezoid rule, or Simpson's rule, to 
evaluate the definite integral in formula (61). Alternatively, any one of the 
functions ~(/o), e (R) or n(R) might be approximated by straight-line segments 
or by parabolic arcs (the approximation of n(R) by straight-line segments would 
amount to assuming a spherically stratified atmosphere). In any case, the re- 
sulting expression for 8^ would have the form of a sum of easily evaluated 
terms. If the number of terms in such a formula is too large for easy calcula- 
tion, it might be feasible to approximate the refractivity profile (or ihe e -func- 
tion) by expressions which are more complex than straight-line segments or 
parabolic arcs, but which are still simple enough so that the resulting integrals 
can be evaluated. In this way a formula might be developed which has a smaller 
number of terms. 

A word of caution is in order here. Before replacing the refractivity profile 
by simpler functions, one should make sure that the corresponding approximate 
e -function is non-negative and satisfies condition (25) . It is not hard to see that 
the various integrals will exist for any piece-wise continuous e -function which 
has these properties. (For example, the proof which was given for the existence 
of the integral in formula (61) applies to any e -function which has these two 
properties, even if it is not continuous.) It follows that these integrals can be 
approximated by using such an approximate e -function, even though the laws of 
geometric optics, and therefore our entire derivation, do not strictly hold (for 
radio waves in the ionosphere) if the refractivity profile is discontinuous. 


Quantities Related to the Arrival Angle 


There are two quantities related to the arrival angle which are widely used; 
the arrival angle at the satellite and the total bending angle. 

The Arrival Angle at the Satellite — This angle can easily be computed by 
using the spherical form of Snell’s law, as given by Equation (13). To carry out 
the calculation we refer to Figure 8 , which should be compared with Figure 1. 

In Figure 8 the vector is a unit vector tangent to the ray path at the satellite 
and all angles are taken positive in the directions indicated. We want to compute 
the angle <>5 between the slant path and the vector s, . In general it seems that 
may be either positive or negative, but in either case we have = 7 + . 
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Figure 8— Diagram for the Calculation of 


Recall also that, for any refractivity profile in the class 8 , < tt /2 (see the 
section entitled "Some General Properties of Ray Paths"); this fact is essential 
for the followii^ derivation. 

From equation (13) it follows that 

Rj sin 0 ^ = cos (a^ + 5^), 
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where Is the value of the index of refraction n at the satellite. If we solve 
this equation for (= -7) we obtain 


5 


sin 


-1 


, rTT" 

L s a 


COS (a^ + 5 ) 


- y . 


where the arc sine, which is the angle 6 ^ , is taken in the first quadrant. We can 
find the angle 7 by using the law of sines; we see that 


s in 7 _ 


or 


7 = sin * 



where the arc sine is again taken in the first quadrant, since it is clear that 
0 <7 < 77/2. We thus obtain the following formula for 5 ^ : 


8 = sin 


-1 




• cos (a^ + 8 ) 


- sin 


-1 
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where both arc sines are taken in the first quadrant. Since 8^ depends on the 
satellite position and on the ref r activity profile, we see that formula (62) gives 
8^ in terms of those same quantities. Note that 8^ - OasR^-'Q^. 

The Total Bending Angle — The total bending angle can be defined as the 
angle through which the vector must be rotated to coincide with the vector 
s^ , with the angle taken positive in the clockwise direction. If we denote the 
total bending angle by , it is easy to see from Figure 8 that 

s g 

= [77 - 7 - (77/2 +0.^)3 + (7 + S,) - iTT /2 - - 8g) , 
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or 


/3 - 8 +8 

'^s g s 


( 63 ) 


Note that /3 - S as R - o^ . 

'^s g s 

A Formula for the Total Bending Angle in Terms of Satellite Altitude and 
the Arrival Angle - Equations (62) and (63) give and /3^ in terms of (= R^ 
plus the satellite altitude), the arrival angle (= eleva- 

tion angle . it is of some interest that we can derive a formula for 8^ which 
depends only on R^ and . if could be accurately measured such a formula 
would be very useful, since /3^ for very distant satellites. Unfortunately 
angular quantities cannot be measured very precisely by radar equipment, so the 
formula would be useful mainly for visual observations. In fact the formula we 
will derive will be specialized (in the next section) to the classical astronomical 
formula for tropospheric bending. In contrast to the case for ,5^ , it does not seem 
possible to obtain formulas for Sg or B^ which depend only on R^ and . 

To derive the formula we first define the partial bending angle /3 at a given 
point along the ray path to be the angle through which the vector must be ro- 
tated to coincide with the vector s at the given point (see Equation (11) for the 
definition of s ) , with the angle taken positive in the clockwise direction. If the 
given point is at the satellite position, then of course ,6 is just /3 ^ . From Figure 
1 we see that 


(- 4 ; + 8 ) ^ 6 ^. 


From our analysis of the ray paths we know that R is a strictly monotonic ally 
increasing function of (//, for a given ray path. It follows that the angles /S , 0 
and 8 are (single-valued) functions of R, for a given ray path. We know also 
that the ray path is differentiable; this can be seen from Equation (20), if we 
recall that any e -function in the class 8 is differentiable and therefore continu- 
ous. It follows that /3 , 0 and 8 are differentiable functions of R. We may there- 
fore differentiate the above equation with respect to R: 


dR " dR dR ’ 
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We can obtain an expression for dt?/dR by differentiating Equation (13): 


dR 


n 

K 


cos a 

a 



n2 r2 


1 



/R \2 

/ e i 2 

— • COS a 

\ R / 


If we insert this into the above equation for d/3/dR together with the expression 
for d0/dR given by Equation (19) some terms cancel out and we obtain the 
simple formula 


d/3 


dR 



- n R cos a 

g e a 


dn 

dR 



(In using Equation (19) recall that 1 + e = formula gives d/3/dR ex- 

plicitly as a function of R, since we assume that the function n(R) is known. If 
we integrate it from R^ to R^ we obtain our formula for ; 




n R cos a 

g e a 


f 


^dR 

dR 


e 


n-'R 


® ^ ' -1 ] cos^ a 

n / \R / 


1- ~ 


-|l/2 
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This formula depends explicitly on R^ and a^, but not on . It is valid for 
any profile n(R) \viiose corresponding c -function is in the class S , for all 
R > R and for 0 < <v/2> For distant satellites the formula can be slightly 

simplified by noting that dn/dR ^ 0 above the sensible atmosphere. If we denote 
the value of R at the altitude where the atmosphere may be considered to end by 
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Rf , then we may w;ite an/dR = 0 for R > R^. If the satellite is extremely dis- 
tant, ~ /3^ and we may then write 


S ~ /3 = -n R cos a 

g * g e a 




1/2 


(65) 


This formula depends explicitly only on the arrival ar^le . It is valid for all 
and for all sufficiently large values of R^. If could be accurately measured 
tlds formula would be a very useful and accurate approximation for S^, for large 
satellite distances. (We could replace by and use the resultii^ expression 
as an approximation for but it probably would not be very accurate except for 
large elevation angles.) Since cannot be measured accurately for radar waves, 
however, formula (65) is useful mainly for visual observations. In this case it 
can be written in a modified form. 

The Classical Formula for the Tropospheric Bending Angle — For visible 
light the ionospheric bending is completely negligible, because of the inverse 
dependence of the refractivity on frequency (see Equation (9)), and we need only 
consider the effect of the troposphere. We may assume that the tropospheric 
refractivity decreases monotonically with altitude, if we disregard local effects 
such as temperature inversions. In other words, we may assume that n(R) is a 
strictly monotonically decreasing function of R throughout the troposphere (see 
Equation (2) and the subsequent discussion). (It follows that R is a (single- 
valued) function of n, and that R(n) also decreases strictly monotonically.) We 
may therefore use n as the variable of integration in Equation (65); the equation 
then takes the form 


/3 = R cos a 

'^s g e a 


f 


dn 



“ 

/n 

/R 1 

1/2 

n2 R 

1 - 


[ —1 ) cos^ a 




\n / 

\R J ® J 
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We have used the facts that n(R^) = n^ and n(R^) ^ 1, where Rj. is taken to be 
the top of the troposphere. Equation (66) is the classical astronomical formula 
for the atmospheric bending of visible light. It is valid for 0 < < tt/ 2. For 

astronomical objects R^ is of course extremely large, so that (66) is a very 
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accurate formula for the elevation angle error A simple derivation of this 
formula is given in Reference 3 (pages 62-65). Note that an accurate measure- 
ment of is all that is needed to compute assuming that the refractivity 
profile is known. In practice, of course, the integral is approximated to obtain 
sim^-le formulas valid for objects not too close to the horizon; for objects near 
the horizon, empirical tables are used (see Reference 3 for a discussion of 
these approximations). 
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APPENDIX 


THE SPERRY REFRACTIVITY PROFILE 


The refractivity profile described in this appendix is referred to in this 
report as the "Sperry profile" because it was used in a report written by the 
Sperry Gyroscope Co. The profile assumes an exponentially decreasing tropo- 
spheric refractivity and a nearly parabolic shape for the F^ layer of the iono- 
sphere. The D and E layers of the ionosphere are neglected. The profile and 
the € -function corresponding to it are merely piecewise continuous » but the 
e-hmction satisfies condition (16) and can therefore be used in the various 
integrals with assurance that they will exist. 

Recall that the refractivity N is defined by 


N = (n - 1) X 10« , 


where n is the index of refraction. The tropospheric profile is assumed to be 
exponentially decreasing: 


N = Nq 


^-h/H 


(Al) 


where N^ (the refractivity at the earth's surface) = 313, h = height (in kilometers) 
above the earth's surface and H (scale height) - 7 kilometers. 

The ionospheric (phase) refractivity is assumed to be given by tiie formula 


N = - 4.03 X 10- * iL , 


(A2) 


where f is the signal frequency in megahertz and M is the nun^r of elrotrons 
per cubic meter. Except for the xinits used for f. Equation (A2) is identical with 
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Equation (9). Equation (A2) is considered to be valid for f> fcritic ii “ ^ 

10 megahertz, where is the maximum electron density in the 
layer, in electrons per cubic meter. If we take 10* ^ electrons/meter ^ as the 
largest likely value for we obtain 8.97 megahertz for fcrituai » which is 
well below any radar or even VHF frequencies. 

The electron density M is assumed to be given by 



M - [1 - (1 - 7 ) 2 ] for 0 < / i 1 . 

(A3) 

and by 




M - sech [7 '4 • (7 - I)] for 7 > 1 , 

(A4) 

where 




h - \ (R - Re) - ^ 

h hrt h - h- 

0 m 0 


ho = height of the base of the layer and h^ = height of the maximum electron 
density in the F^ layer. Note that 


scch X - 

cosh X 


2 e- 


e* + e 


+ e 


- 2x 


the next to the last form is probably the most convenient for numerical calcu- 
lation, while the last form clearly shows the decreasing exponential behavior 
of seek X for large x. 

The Sperry profile can now be precisely given as follows: 

• For 0 < h < 40 kilometers, (Al) holds; 

• For 40 km 2 hs h(,, N = 0; 

• For hg < h< h^ , (A3) holds; 

• For h^ < h < 2,000 km, (A4) holds; 

• For h > 2,000 km, N = 0. 
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Note that the ionospheric (i.e„ Fj -layer) N-profile is parabolic below h^ ; above 
h„ it is nearly parabolic, but trails off to an exponential shape at great heights. 
Note also that the tropospheric profile is fixed, while the ionospheric profile de- 
pends (for a given frequency) on three parameters: , h^ and . 

If, instead of N, we work with e , which we recall is defined by the relation 


€ 


n 1 + N. 

- 1 = 2 . 

" 1 + 10‘* N 


- 1 , 


it is clear that the € -function corresponding to the Sperry N-profile is defined 
as follows: 

• For < R < (R^ + 40) kilometers. 


e (R) = 


1 + 10’® Nq 


- 1 ; 


• For (R^ + 40) km < R < R^ + 

^ (R) = 10’® % ; 


• For R, +ho < R <R^ + , 


e (R) = 


1 + 10'® Njj 

1 - 4. 03 X 10"“ [1 - ( 1 - y)2] 

f2 


• For R +h<R< (R + 2,000) km, 

cm 


e (R) = 


1 + 10 '® \ 

1 - 4.03 X 10'^* sech [tt ''4 * ( / - 1)] 
f2 


- l! 
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• For R > (R^ + 2,000) km. 


e(R) = 10-« Nq . 

The corresponding Rinction e (p) is easily obtained by using the transformation 
P= R, /R. We note that the function e (R) is piecewise continuously differen- 
tiable, but that it is discontinuous at R = R^ +40 and R = R^ + 2,000. In other 
words, it is discontinuous at the top of the troposphere and at the top of the 
ionosphere. These discontinuities cause no dlMculty, however, because e (R) 
is clearly always positive (except that ) = 0), and we will show that it satis- 
fies condition (25). It is not hard to see that there must therefore exist functions 
arbitrarily close toe (R), but continuously differentiable everywhere, so that they 
are in the class 8 ^ • Our entire analysis is valid for any of these fuiutions, and 
the Sperry 6 -function may be regarded as a simple approximation to any of them. 
All of the definite integrals in the various equations and formulas are therefore 
guaranteed to exist if the Sperry profile is used. 

It remains only to show that e (R) satisfies condition (25). To show that 
€(R) is bounded is easy; we just note tiiat the maximum value of e (R) occurs 
at R = R ^ + h^ , at which point y = 1, and therefore 

1 + 10-« Nq 

6(R) < 1 

1 - 4.03 X 10-^1 J!!l 
f2 


for ail R. Further, if we assume that 100 megahertz is the smallest frequency 
likely to be used (corresponding to VHF frequencies), and if we assume 10^^ 
electrons/M^ as the largest value of , we get the bound 




, , 1 + 0.000313 
^ 1 - 0.00403 


= 0.00436, 


which holds for any realistic profile and for any signal frequency in the VHF 
band or higher. 

To show that e (R) < (1 -cr ) (R - R^ )/l^ for someo-8 (0,1) is more trouble- 
some, but still not di^icult. We must first look at the tropospheric profile. For 
R^l R< (R^ + 40) km we have 
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e (R) 


1 + a 
1 + aE 


- 1 , 


and 


d e _ a(l + a)E 
U(l+aE')2 


dR^ H2 (1+ aE)3 


where we have written E = exp [ -(R - Re )/H] and a - 10 ® No for brevity. 
Since E ^ 1 for R - Re , we have 1-aE^ l-a=l-10 313 > 0; thus the 

second derivative is negative, and of course the first derivative is positive. It 
follows that the maximum slope of ^ (R) occurs at R = R^ , so that we may 
write 


f ( 


R) < ^ 

dR 


• (R-Re) 


aR 

e 

H(l+a) 


(R - Re) 

P 


(R - Re) 

^ (1 - 0.715) • 


-L(R) , 


where we have used the value = 6,378 km. This shows that the required con- 
dition is satisfied by the tropospheric profile, with the value cr = 0.715. 

Next we look at the region between the top of the troposphere and the bottom 
of the ionosphere. We find that 


L(R^ + 40) = 0.00179 > e (R^ + 40) = 0.000313; 


thus the e -function is bounded above by L(R) throughout this region, and the 
same is obviously also true of the region above the top of the ionosphere. 
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Finally we will show that t (R) < L(R) throughout the part of the ionosphere 
below the Fj maximum; the same will then be true of the rest of the ionosphere, 
because £ (R) decreases above the F^ maximum. For +h^^< R< R^ + h^ 
we find that 


£ (R) 


1 - b[l - (1 - 7 ) 2 ] 


d£ _ 2b(l + a)(l - 7) 

” {1 - b[l - (1 = 

d^e 2b(l + a)[-l+b + 3b(l - 7)2] 
d7^ {1 - b[l - (1 - 7)2] }3 


where we have written a = 10~® Nq and b = 4.03 x 10““ /f2 for brevity. If 

we again take 10^2 electrons/M^ and 100 megahertz as limiting values for 
and f, respectively, we find that the largest possible value of b is given by b = 
4.03 X 10“*^ X 10 ^ 2/104 = 0.00403. It follows that [-1 + b + 3b (1 - 7 ) 2 ] <- 1 + 

second derivative is n^ative. Of more interest for our 
present purpose is the fact that the first derivative is positive. (Recall that 
0 < 7 ^ 1 in the portion of the ionosphere below the F^ maximum). It follows that 
it would be sufficient (but not necessary) to show that £ (R^ + h^) < L (R + h^). 
This condition will be satisfied if (1 +a)/(l - b^^^ ) ~ 1 < 0.285 h^^/R (recall 
that 7=1 for R = R^ + h^); i.e., if h^ > 98 kilometers. Since the layer is 
usually considered to begin at the considerably higher altitude of about 140 
kilometers (see for example Reference 9), we may safely conclude that the 
linear function L(R) is in fact a majorant for the Sperry profile, for any real- 
istic values of the parameters h^, h^ and • 

In the computer program, which will be described in a companion report, 
there is a subroutine which numerically tests the e -function to determine whether 
it satisfies condition (25). The Sperry profile was tested by this routine and 
was found to satisfy the condition for arbitrarily chosen, but realistic, values of 
the parameters hg , h^ and . Several sets of values for these parameters 
were used in the computer runs, but representative values, used for most of the 
runs, are: h„ = 240 km, h = 300 km and M = 5.2 x 10“ electrons per 
cubic meter. 
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